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L Introduction 



In coding theory, the so-called linear programming method, introduced by Phi- 
lippe Delsarte in the seventies |16| as proved to be a very powerful method to 
solve extremal problems. It was initially developed in the framework of association 
schemes and then extended to the family of 2-point homogeneous spaces, including 
the compact real manifolds having this property (see ifTSl . 1241 . lfT3l Chapter 9]). 
Let us recall that a 2-point homogeneous space is a metric space on which a group 
G acts transitively, leaving the distance d invariant, and such that, for {x,y) € X^, 
there exists g € G such that {gx, gy) = {x', y') if and only if d{x, y) = d{x', y'). 
The Hamming space Hn and the unit sphere of the Euclidean space 5""^^ are core 
examples of such spaces which play a major role in coding theory. To such a 
space is associated a sequence of orthogonal polynomials {Pk)k>o such that, for 



These inequalities can be understood as linear constraints on the distance distribu- 
tion of a code and are at the heart of the LP method. 

The applications of this method to the study of codes and designs are numerous: 
very good upper bounds for the number of elements of a code with given minimal 
distance can be obtained with this method, including a number of cases where this 
upper bound is tight and leads to a proof of optimality and uniqueness of certain 
codes, as well as to the best known asymptotic bounds (see |[T6l . 1301 . ll24l . lfT3l 
Chapter 9], [28|). 

In recent years, the development of the theory of error correcting codes has in- 
troduced many other spaces with interesting applications. To cite a few, codes 
over various alphabets associated to various weights, quantum codes, codes for the 
multi antenna systems of communications involving more complicated manifolds 
like the Grassmann spaces, have successively focused attention. For these spaces 
there was a need for a generalization of the classical framework of the linear pro- 
gramming method. This generalization was developed for some of these spaces, 
see 1441 . B31 . El . |[37ll . It turns out that in each of these cases, a certain sequence 
of orthogonal polynomials enters into play but unlike the classical cases, these 
polynomials are multivariate. 

Another step was taken when A. Schrijver in BOl succeeded to improve the 
classical LP bounds for binary codes with the help of semidefinite programming. 
To that end he exploited SDP constraints on triples of points rather than on pairs, 
arising from the analysis of the Terwilliger algebra of the Hamming scheme. His 
method was then adapted to the unit sphere IH in the framework of the represen- 
tations of the orthogonal group. The heart of the method is to evidence matrices 



all C CX, 



Pk{d{c,c'))>0. 
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Zfc(x, y, z) such that for all C C X, 

Zkic,c',c")hO. 

{c,c',c")£C3 

Another motivation for the study of SDP constraints on /c-tuples of points can 
be found in coding theory. It appears that not only functions on pairs of points such 
as a distance function d{x, y) are of interest, but also functions on A;-tuples have 
relevant meaning, e.g. in connection with the notion of list decoding. 

In these lecture notes we want to develop a general framework based on har- 
monic analysis of compact groups for these methods. In view of the effective ap- 
plications to coding theory, we give detailed computations in many cases. Special 
attention will be paid to the cases of the Hamming space and of the unit sphere. 

Section 2 develops the basic tools needed in the theory of representations of fi- 
nite groups, section 3 is concerned with the representations of compact groups and 
Peter Weyl theorem. Section 4 discusses the needed notions of harmonic analy- 
sis: the zonal matrices are introduced and the invariant positive definite functions 
are characterized with Bochner theorem. Section 5 is devoted to exphcit computa- 
tions of the zonal matrices. Section 6 shows how the determination of the invariant 
positive definite functions leads to an upper bound for codes with given minimal 
distance. Section 7 explains the connection with the so-called Lovasz theta num- 
ber. Section 8 shows how SDP bounds can be used to strengthen the classical LP 
bounds, with the example of the Hamming space. 

1.1. Notations: for a matrix A with complex coefficients. A* stands for the trans- 
posed conjugate matrix. A squared matrix is said to be hermitian if A* = A and 
positive semidefinite if it is hermitian and all its eigenvalues are non negative. This 
property is denoted A ^ 0. We follow standard notations for sets of matrices: the 
set of n X m matrices with coefficients in a field K is denoted K"'^'^; the group 
of n X n invertible matrices by Gl(i^"); the group C/(C") of unitary matrices, 
respectively 0{W^) of orthogonal matrices is the set of matrices A £ G1(C"), re- 
spectively A e G1(]R") such that A* = A'^. The space C"'''" is endowed with 
the standard inner product {A,B) = Trace{AB*) = J2i j The number 

of elements of a finite set X is denoted card(X) of 

2. Linear representations of finite groups 

In this section we shortly review the basic notions of group representation theory 
that will be needed later. There ai^e many good references for this theory e.g. 141.1 . 
or Il38l which is mainly devoted to the symmetric group. 

2.1. Definitions. Let G be a finite group. A (complex linear) representation of G 
is a finite dimensional complex vector space V together with a homomorphism p: 

p:G^ Gl{V) 

where G1(F) is the general linear group of V, i.e. the set of linear invertible trans- 
formations of V. The degree of the representation {p, V) is by definition equal to 
the dimension of V. 

Two representations of G say (p, V) and {p', V') are said to be equivalent or 
isomorphic if there exists and isomorphism u : V ^ V such that, for all g £ G, 

p'{g) = up{g)u~^. 
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For example, the choice of a basis of V leads to a representation equivalent to 
{p, V) given by {p', C^) where d = diin(F) and p'{g) is the matrix of p{g) in the 
chosen basis. In general, a representation of G such that V = is called a matrix 
representation. 

The notion of a G-module is equivalent to the above notion of representation and 
turns out to be very convenient. A G-module, or a G-space, is a finite dimensional 
complex vector space V such that for all g ^ G, v ^ V , gv £ V is well defined 
and satisfies the obvious properties: Iv = v, g{hv) = {gh)v, g{v + w) = gv + gw, 
g{Xv) = \{gv) for g,h e G, v,w e V, X e C In other words, V is endowed 
with a structure of C[G] -module. One goes from one notion to the other by the 
identification gv = p{g){v). The notion of equivalent representations corresponds 
to the notion of isomorphic G-modules, an isomorphism of G-modules being an 
isomorphism of vector spaces u : V ^ V' such that u{gv) = gu{v). Note that 
here the operations of G on F and V' are denoted aUke, which may cause some 
confusion. 

2.2. Examples. 

• The trivial representation 1: = C and gv = v. 

• Permutation representations: let X be a finite set on which G acts (on the 
left). Let Vx '■= (BxexCcx- A natural action of G on Vx is given by 
g^x = Ggx, and defines a representation of G, of degree The matrices 
of this representation (in the basis {e^;}) are permutation matrices. 

- The symmetric group Sn acts on X = {1,2,..., n}. This action 
defines a representation of degree n of S^. 

- For allw,l < w < n, Sn acts on the set Xyj of subsets of {1, 2, . . . , n} 
of cardinal w. In coding theory an element of Xyj is more Ukely 
viewed as a binary word of length n and Hamming weight w. The 
spaces Xyj are called the Johnson spaces and denoted J^. 

• The regular representation is obtained with the special case X = G with 
the action of G by left multiplication. In the case G = Sn it has degree 
n\.. It turns out that the regular representation contains all building blocks 
of all representations of G. 

• Permutation representations again: if G acts transitively on X, this action 
can be identified with the left action of G on the left cosets G/H = {gH : 
g & G} where H = Stab(xo) is the stabihzer of a base point. 

- The symmetric group Sn acts transitively on X = {1,2, ... ,n} and 
the stabilizer of one point (say n) can be identified with the symmetric 
group Sn-i acting on {1, . . . , n — 1}. 

- The action of Sn on is also transitive and the stabilizer of one point 
(say i"'0"~"') is the subgroup S'ji x S'j^+i „} isomorphic to 

Sw ^ Sn^w 

- The Hamming space = {0, 1}" = F2 affords the transitive action 
of G = T X Sn where T is the group of translations T = {tu ■ u & 
Hn}, tu{v) = u + v and Sn permutes the coordinates. The stabihzer 
of 0" is the group of permutations 5„. 

• Another way to see the permutation representations is the following: let 



C{X) -{f-.X^C} 
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be the space of functions from X toC. The action of G on X extends to a 
structure of G-module on C{X) given by: 

For the Dirac functions 5y {Sy{x) = 1 if x = y, otherwise), the ac- 
tion of G is given by gSy = Sgy thus this representation is isomorphic to 
the permutation representation defined by X. This apparently more com- 
plicated presentation of permutation representations has the advantage to 
allow generalization to infinite groups acting on infinite spaces as we shall 
encounter later. 

2.3. Irreducibility. Let F be a G-module (respectively a representation (p, V) of 
G). A subspace W C V is said to be G-invariant (or G-stable, or a G-submodule, 
or a subrepresentation of {p, V)), H gw e W (respectively p{g){w) G W) for all 
g e G,w gW . 

Example: V = Vq and W = Cec with ec = Z^^gg ^9- restriction of the 
action of G to is the trivial representation. 

A G-module V is said to be irreducible if it does not contain any subspace W, 
W 7^ {0}, V, invariant under G. Otherwise it is called reducible. 
Example: The G-invariant subspaces of dimension 1 are necessarily irreducible. 
If G is abelian, a G-module of dimension greater than 1 cannot be irreducible, 
because endomorphisms that pairwise commute afford a common basis of eigen- 
vectors. 

The main result is then the decomposition of a G-module into the direct sum of 
irreducible submodules: 

Theorem 2.1 (Maschke's theorem). Any G-module V ^ {0} is the direct sum of 
irreducible G -submodules Wi, . . . , W^: 

(1) V = Wi®W2®---®Wk. 

Proof. By induction, it is enough to prove that any G-submodule WofV affords 
a supplementary subspace which is also G-invariant. The main idea is to construct 
a G-invariant inner product and then prove that the orthogonal of W for this inner 
product makes the job. 

We start with an inner product {x, y) defined on V. There are plenty of them 
since F is a finite dimensional complex vector space. For example we can choose 
an arbitrary basis of V and declare it to be orthonormal. Then we average this inner 
product on G, defining: 

{x,yy := ^{gx,gy). 

It is not difficult to check that we have defined a inner product which is G-invariant. 
It is also easy to see that 

:= {veV : {v, w)' = for all weW] 

is G-invariant, thus we have the decomposition of G-modules: 

V = W®W^ 



□ 
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It is worth to notice that the above decomposition may not be unique. It is clear 
if one thinks of the extreme case G = {1} for which the irreducible subspaces are 
simply the one dimensional subspaces of V. The decomposition of V into the direct 
sum of subspaces of dimension 1 is certainly not unique (if dim(F) > 1 of course). 
But uniqueness is fully satisfied by the decomposition into isotypic subspaces. In 
order to define them we take the following notation: let 7^ be a complete set of 
pairwise non isomorphic irreducible representations of G. We have seen that any 
G-module affords a G-invariant inner product so the action of G on ii is expressed 
by unitary matrices in a given orthonormal matrix of R. According to the context 
we view R either as a G-module or as a homomorphism g i— )• R{g) G U (C"). It 
will turn out that there is only a finite number of them but we have not proved it yet. 
The isotypic subspace Ir of V associated to i? G 7^ is defined, with the notations 
of O.by: 

(2) Ir := ew^c^RWi. 

Theorem 2.2. Let R & TZ. The isotypic spaces Ir do not depend on the decom- 
position ofV as the direct sum of G -irreducible subspaces. We have the canonical 
decomposition 

V = ®RenlR. 

Any G-subspace W C V such that W R is contained in Ir and any G- 
irreducible subspace o/Ir is isomorphic to R. A decomposition into irreducible 
subspaces o/Ir has the form 

iR = Wi(B---e 

with Wi ~ R. Such a decomposition is not unique in general but the number uir 
does not depend on the decomposition and is called the multiplicity of R in V. 

Moreover, if V is endowed with a G-invariant inner product, then the isotypic 
spaces are pairwise orthogonal. 

Proof. We start with a lemma which points out a very important property of irre- 
ducible G-modules. 

Lemma 2.3 (Schur Lemma). Let Ri and R2 two irreducible G-modules and let 
(p : Ri ^ R2 be a G -homomorphism. Then either ip = or (p is an isomorphism 
of G-modules. 

Proof. The subspaces ker ip and im p are G-submodules of respectively Ri and 
i?2 thus they are equal to either {0} or Ri. □ 

We go back to the proof of the theorem. We start with the decomposition ([T]) 
of V and the definition Q of Ir, a priori depending on the decomposition. Let 
W C V,a. G-submodule isomorphic to R. We apply Lemma [231 to the projections 
pw^ and conclude that either PWi{W) = {0} or PWi{W) = Wi and this last case 
can only happen if W Wi. It proves that W C Ir and that a G-irreducible 
subspace of Ir can only be isomorphic to R. It also proves that 

Ir= E ^ 

WcV,W~R 

hence giving a characterization of Ir independent of the initial decomposition. 
The number niR must satisfy dim(XR) = uir dim(i?) so it is independent of the 
decomposition of Ir. 
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If V is equipped with a G-invariant inner product, we consider orthogonal pro- 
jections. Schur Lemma shows that Pw{W') = {0} or = W if W and W are 
irreducible. Thus if they are not G-isomorphic, W and W must be orthogonal. 

□ 

2.4. The algebra of G-endomorphisms. Let be a G-module. The set of G- 
endomorphisms of V is an algebra (for the laws of addition and composition) de- 
noted EndciV). The next theorem describes the structure of this algebra. 

Theorem 2.4. //F ~ e/^e7^-R"'^, then 

EndG(l^) ^ n C™«^™«. 
Ren 

Proof. The proof is in three steps: we shall assume first V = Ris irreducible, then 
V ~ i?*", then the general case. Schur Lemma [231 is the main tool here. 

If V is irreducible, let ip G Endc{V). Since F is a complex vector space, if has 
got an eigenvalue A. Then — A Id is a G-endomorphism with a non trivial kernel 
so from Schur Lemma (/? — A Id = 0. We have proved that 

EndG(^) = {AId,A G C} ~ C. 

We assume now that V ~ ii™ and we fix a decomposition V = Wi © • • • © Wm. 
For all 1 < i < j < m, let Uj^i : Wi — )• Wj an isomorphism of G-modules such 
that the relations 

Ukj ° Uj^i = Uk,i and Ui^i = Id 

hold for all i, j, k. Let Lp G EndG(^); we associate to Lp an element of C™^™ in 
the following way. From previous discussion of the irreducible case it follows that 
for all i,j there exists Ojj G C such that, for all v G Wi, 

PWj ° ^(v) = aj^iUj^i{v). 

The matrix A = {aij) is the matrix associated to ip. The proof that the mapping 
9? I— )• j4 is an isomorphism of algebras carries without difficulties and is left to the 
reader. 

In the general case, V = ®Ren'^R- Let (p G EndG(V^)- It is clear that p{1r) C 
Ir thus 

EndG(F) = ©Rg7eEndG(XR) 
and we are done. □ 

It is worth to notice that EndG(^) is a commutative algebra if and only if all the 
multiplicities m/j are equal to either or 1. In this case we say that V is multiplicity 
free. It is also the unique case when the decomposition into irreducible subspaces 
([T]) is unique. 

2.5. Orthogonality relations. Another important result which is a consequence 
of Schur lemma is the orthogonality relations between the matrix coefficients of 
the elements of TZ: 

Theorem 2.5. For R e TZ, let dn := dim(fi). For all R,S eU, i, j, k, I, 

{Ri,j,Sk,l) = -^SR,S^i,kSj,l- 
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Proof. For A G C'^^^'^s^ let 

' ' sec 

This matrix satisfies R{g)A' = A'S{g) for all g £ G. In other words it defines 
an homomorphism of G-modules from (C"'^, 5) to (C"'^, R). Schur lemma shows 
that if S ^ R, A' = and if S = R, A' = Aid. Computing the trace of A' 
shows that A = Trace{A) / cIr. Taking A = Eij the elementary matrices, with the 
property that S{g)^^ = S{g)*, leads to the announced formula. □ 

2.6. Characters. The character of a representation (p, V) of G is the function 
Xp ■ G ^ C defined by 

Xpig) = Trace(p(ff)). 
As a consequence of the standard property of traces of matrices TTace{AB) = 
Trace (i? A), the character of a representation only depends on its equivalence class, 
and it is a complex valued function on G which is constant on the conjugacy classes 
of G (such a function is called a class function). The inner product of any two x, 
tp gC{G) is defined by 



1 



gee 

We have the very important orthogonality relations between characters: 

Theorem 2.6 (Orthogonality relations of the first kind). Let x '^nd x' be respec- 
tively the characters of two irreducible representations {p,V) and {p',V') of G. 
Then 

[l ifp^p' 
I otherwise. 



{x,x') 



Proof. It is a straight forward consequence of Theorem I2.5[ since the trace of a 
representation is the sum of the diagonal elements of any equivalent matrix repre- 
sentation. □ 

A straightforward consequence of the above theorem is that {xp, Xr) = 
for all i? € TZ. This property is a very convenient tool to study the irreducible 
decomposition of a given representation {p, V) of G; in particular it shows that a 
representation is irreducible if and only if its character x satisfies {x,x) = 1- Iri 
the case of the regular representation it leads to the following very important result: 

Theorem 2.7. [Decomposition of the regular representation] 

Proof. Compute the character of the regular representation. □ 

A consequence of the above theorem is the finiteness of the number of irre- 
ducible representations of a given finite group, together with the formula 

|G|= 5](dim(i?))2 
rgtz 

which shows e.g. completeness of a given set of irreducible G-modules. 

A second consequence of the orthogonality relations is that a representation of 
G is uniquely characterized up to isomorphism by its character. 
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Theorem 2.8. 

{p,V)^{p',V')^Xp = Xp'- 

Proof. If xp = Xp' ' the multiplicities of an irreducible representation of G are the 
same in V and V' , hence V ~g V . □ 

Let us denote by R{G) the subspace of elements of C{G) which are constant on 
the conjugacy classes Ci, . . . , of G. The dimension of R{G) is obviously the 
number s of conjugacy classes of G. We have seen that the characters XR of the 
irreducible representations of G belong to R{G) and form an orthonormal family. 
It turns out that they in fact form a basis of R{G), which in other words means 
that the number of irreducible representations of G is exactly equal to its number 
of conjugacy classes. 

Theorem 2.9. 

R{G) = (BB.GnCxR. 

Proof. It is clear that C(G) = C[G]5e. Thus EndG(C(G)) ~ C[G]. In particular, 
the center of Endc{C{G)) is isomorphic to the center Z{C[G]) of C[G]. It is easy 
to verify that the center of C[G] is the vector space spanned by the elements Aj := 
SgGC 9 associated to each conjugacy class Cj of G, thus Z{C[G]) is of dimension 
s the number of conjugacy classes of G. On the other hand, as a consequence of 
Theorem [277] and Theorem f2A[ we have EndG(C(G)) ~ 11^67?^'^'*'"^'' where 
dn = dim(i2). Thus the center of EndG(C(G)) is isomorphic to C'^' and we 
have proved that the number of G-irreducible modules is equal to the number of 
conjugacy classes of G. □ 

Remark 2.10. There is not in general a natural bijection between the set of con- 
jugacy classes of G and the set of its irreducible representations. However, in the 
special case of the symmetric group Sn, such a correspondance exists. The conju- 
gacy classes are naturally indexed by the partitions Xofn and to every partition X 
ofn is associated an irreducible module S''^ also called a Specht module (see \3Sl ). 

2.7. Induced representation and Frobenius reciprocity. Induction is a way to 
construct representations of a group G from representations of its subgroups. Look- 
ing at the irreducible subspaces of representations that are induced from subgroups 
is a very convenient way to find new irreducible representations of a group G. In- 
duction is an operation on representations which is dual to the easier to understand 
restriction. If F is a G-module and H h a subgroup of G, the restriction Res^(T^) 
is simply the space V considered as a C[i7] -module. If V is an //-module, we 
define Indg(y) to be the €[0] -module 

Ind|(y) := C[G] (^ciH] V. 

Here we exploit the bi-module structure of C[G] (the tensor product over C[H] 
means that Xp (>i) v = X (>S> pv when p G C [//]). A more explicit (but less intrisic) 
formulation for this construction is the following: let {xi, . . . , xt} be a complete 
system of representatives ofG/H, so that G = xiH VJ ■ ■ ■ VJ xtH. Then 

Indg(y) = eU^XiV 

where the left action of G is as follows: for all i, there is j and h £ H both 
depending on g such that gxi = Xjh. Then gxiV := Xj{hv) where hv ^ V. K 
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third construction of Ind^(y) is the following: 



Indg(F) = {f:G^V such that f{gh) = h-'f{g)}. 

The equivalence of these three formulations is a recommended exercise ! 

Example: The permutation representation of G on X = G/H is nothing else than 
the induction of the trivial representation of H. In short, C(X) = Ind^ 1. 

Since the induction of two isomorphic i7-modules leads to isomorphic G-modules 
and similarly for the restriction, these operations act on the characters thus we de- 
note similarly Res^ X' Ind^ x the characters of the corresponding modules. 

Lemma 2.11. Let x a character of H. The induced character Ind^ x given 
by the formula: 

Proof. We take a decomposition Ind^(y) = xiV (B- ■ -(BxtV where {xi, ... ,xt} 
are representatives of G/H. Since gxiV = xjhv with the notations above, g{xiV) C 
XjV and the block XiV will contribute to the trace of x gx only when j = i, 
which corresponds to the case when x~^gxi € H. Then, the multiplication by g 
on XiV acts like the multiplication by h = x^^gxi on V. Thus we have 

l<i<t 

x~^gx&H 

□ 



The duality between the operations of restriction and induction is expressed in 
the following important theorem: 

Theorem 2.12 (Frobenius reciprocity). Let H be a subgroup ofG and let x (^nd ijj 
be respectively a character of H and a character ofG. Then 

(Indgx,^') = (x,ResgV')- 
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Proof. Let x : G — )• C be defined by: x{g) = xid) 9 ^ H and x{9) = if 
g ^ H (of course x is not a character of G). We compute (Ind^ Xi 

(Indg X, V') = ^ E I^dg X(5)^ 

' " ' geG xeG 



' " ' a;GG geG 

1 



xGG g'GG 

" ' xgG s'gG 

^5]x(/i)VK^=(x,Resgv). 



□ 



2.8. Examples from coding theory. In coding theory we are mostly interested in 

the decomposition of C{X) under the action of G = Aut(X) for various spaces 
X. We recall that the action of G on / G C{X) is given by {gf){x) = f{g~^x). 
The space C{X) is endowed with the inner product 

which is G-invariant. 



2.8. 1. The binary Hamming space if^: recall that G = TyiSn- Let, for y G Hn, 
Xy & C{Hn) be defined by Xyi^) = (—1)^'^- The set {xy,y £ Hn} is exactly 
the set of irreducible characters of the additive group F2 , and form an orthonormal 
basis of C{Hn). The computation of the action of G on Xy shows that for a G Sn, 
<^Xy = Xa{y) and for tu G T, tuXy = {-'^T'^Xy Let, for < A; < n, 

-ffe ■~-^y,wt{y)=k ^Xy 

Thus Pk is a G-invariant subspace of C{Hn) of dimension (^) and we have the 
decomposition 

(3) C{Hn) = Po±Pl±--- ±Pn. 

The computation {xPkiXPk) = 1 where xPk t^^^ character of the G-module P^ 
shows that these modules are G-irreducible. 

Now we introduce the Krawtchouk polynomials. The element Z^^ := Ylwt{y)=k Xy 
of C{Hn) is S'n-invariant. In other words, Zi^ix) only depends on wt{x). We define 
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the Krawtchouk polynomial for < A; < n by 

(4) Kk{w) : = Zk{x) = ^ {-If-y where wt{x) = w 



wt(y)=k 

w 

(5) ^ 



w \ I n — w 
k — i 



1=0 \ / \ 

We review some properties of these polynomials: 

(1) deg(i^fe) = k 

(2) Kk{0) = il) 

(3) Orthogonality relations: for all < A; < Z < n 



i,t(:)^.(-)A-,c.)=4,,(';) 



The last property is just a reformulation of the orthogonality of the Z). G Pk, since, 
if /) /' € C{Hn) are -invariant, and f{w) := f{x), wt{x) = w, 



xeHn 

= ^E(:)/w/'W. 

w=0 ^ ' 



The above three properties characterize uniquely the Krawtchouk polynomials. 

Let C C Hn be a binary code. Let Ic be the characteristic function of C. The 
obvious inequalities hold: 

(6) 0<fe<n, ^ (lc,X2/)'>0. 
Since the decomposition of Ic over the basis Xy reads 

the above inequalities are indeed reformulations of the non negativity of the squared 
norm of the projections pp^(lc). They express in terms of the Krawtchouk poly- 
nomials: 

(7) 0^^^^' ^ X Kk{dn{x,x'))>^ 

or equivalently in terms of the distance distribution of the code C: if 

A^(C) ■= p|{(x,x') G : dnix.x') = w}\ 

then 

w=0 

These inequalities are the basic inequalities involved in Delsarte linear program- 
ming method. We shall encounter similar inequalities in a very general setting. 
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In the special case when C is linear, we have 

\C\ 

so that we recognise the identity 

wt{y)=k w=0 

to be the Mac Williams identity 

1 " 

The more general g-ary Hamming space affords similar results; it is treated in 
1531 

2.8.2. The Johnson spaces J^: the group is G = Sn- Here, we shall see at work 
a standard way to evidence G-submodules as kernels of G-endomorphisms. For 
details we refer to [17] where the q-Johnson spaces are given a uniform treatment. 
We introduce the applications 

f ^ S{f) : S{f){x) := J2 f^y) 
y&J^, xCy 

and 

^p:c{j:-')^c{j:) 

y&Jn~^, ycx 

Both of these applications commute with the action of G. They satisfy the follow- 
ing properties: (/, tp{f')) = ((^(/), /')> V' is injective and S is surjective. Therefore 
the subspace of C( J^): 

i^u, := ker S 

is a G-submodule of dimension (^) — (J^^) and we have the orthogonal decom- 
position 

By induction we obtain a decomposition 

which can be proved to be the irreducible decomposition of C( J^) (see l5.3.ll) . 
3. Linear representations of compact groups 



In this section we enlarge the discussion to the representation theory of compact 
groups. For this section we refer to lfT2l . 
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3. 1. Finite dimensional representations. The theory of finite dimensional repre- 
sentations of finite groups extends nicely and straightforwardly to compact groups. 
A finite dimensional representation of a compact group G is a continuous homo- 
morphism p : G — )• Gl(y) where y is a complex vector space of finite dimension. 

A compact group G affords a Haar measure, which is a regular left and right 
invariant measure. We assume this measure to be normalized, i.e. the group G has 
measure 1. With this measure the finite sums over elements of a finite group can be 
replaced with integrals; so the crucial construction of a G-invariant inner product 
in the proof of Maschke theorem extends to compact groups with the formula 

:= / {gx,gy)dg. 
Jg 

Hence Maschke theorem remains valid for finite dimensional representations. We 
keep the notation TZ for a set of representatives of the finite dimensional irreducible 
representations of G, chosen to be representations with unitary matrices. A main 
difference with the finite case is that TZ is not finite anymore. 

3.2. Peter Weyl theorem. Infinite dimensional representations immediately oc- 
cur with the generalization of permutation representations. Indeed, if G acts con- 
tinuously on a space X, it is natural to consider the action of G on the space C{X) 
of complex valued continuous functions on X given by {gf ){x) = f{g^^x) to be 
a natural generalization of permutation representations. A typical example of great 
interest in coding theory is the action of G = 0(M") on the unit sphere of the 
Euclidean space: 

:= {x G M" : x-x = l}. 

The regular representation, which is the special case C(G), with the left action of 
G on itself, can be expected to play an important role similar to the finite case. It 
is endowed with the inner product 

(/,/') := / f{9)7i9)dg. 
JG 

For R ^ TZ, the matrix coefficients g — )• Rij{g) belong to unitary matrices. The- 
orem 12.51 establishing the orthogonality relations between the matrix coefficients 
of the elements of TZ remains valid; thus they form an orthogonal system in C(G). 
The celebrated Peter Weyl theorem asserts that these elements span a vector space 
which is dense in C(G) for the topology of uniform convergence. 

Theorem 3.1. [Peter Weyl theorem] The finite linear combinations of the fiinctions 
Rij are dense in C{G) for the topology of uniform convergence. 

Proof. We give a sketch of the proof: 

(1) If y is a finite dimensional subspace ofC{V) which is stable by right trans- 
lation (i.e. by gf{x) = f{xg)) and / G F, then / is a linear combination 
of a finite number of the Ri^: according to previous discussion, there is 
a decomposition V = Wi © • • • © Wn such that Wk is irreducible. If 
Wk — R, there exists a basis ei, . . . , e^^ of Wk in which the action of G 
has matrices R. Explicitly, 

dR 

ej{hg) = ^Ri,j{g)ei{h). 
1=1 
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Taking = 1, we obtain Cj = ej(l)i?ij. 

(2) The idea is to approximate / G C{G) by elements of such subspaces, 
constructed from the eigenspaces of a compact selfadjoint operator. We 
introduce the convolution operators: let cj) G C{G), 

T^{f){9) = icP*f)ig)= [ ct>[gh-^)f[h)dh. 

JG 

(3) Since G is compact, / is uniformly continuous; this property allows to 
choose (f) such that ||/ — T^{f)\\ao is arbitrary small. 

(4) The operator is compact and can be assumed to be selfadjoint. The 
spectral theorem for such operators on Hilbert spaces (here asserts 
that the eigenspaces V\ := {/ : T^f = A/} for A / are finite dimen- 
sional and that the space is the direct Hilbert sum ®a^a- For t > 0, the 
subspaces Vt := |A|>t have finite dimension (i.e. there is only a finite 
number of eigenvalues A with |A| > t > 0). 

(5) The operator commutes with the action of G by right translation thus 
the subspaces Vx are stable under this action. 

(6) Let fx be the projection of / on Va- The finite sums ft := J2\x\>tf>^ 
converge to / — /o for the L^-norm when t — 0. 

(7) Moreover, for all / e C{V), \\T^if)\\oo < ll-Z-llooH/lb. Thus, T^ft) 
converges uniformly to T^{f-fo) = T^{f)- Finally, T^{ft) e Vt and Vt is 
finite dimensional and invariant under the action of G by right translations, 
thus by (1) T^{ft) is a linear combinations of the 

□ 

If dR = dim(i?), the vector space spanned by {Rij,i = l,...,dji} is G- 
invariant and isomorphic to R. So Peter-Weyl theorem means that the decomposi- 
tion of the regular decomposition is 

C{G) =±Ri^n 

where Zr ~ W^^ , generalizing Theorem 12.71 (one has a better understanding of 
this decomposition with the action of G x G on G given by {g,g')h = ghg'~^. For 
this action C{G) = ®/^e7^i^ ® R* where R* is the contragredient representation, 
and i? (g) i?* is G X G-irreducible). 

Since uniform convergence is stronger than convergence, we also have as 
a consequence of Peter Weyl theorem that the matrix coefficients Rij (suitably 
rescaled) form an orthonormal basis of L'^{G) in the sense of Hilbert spaces. 

A slightly more general version of Peter Weyl theorem deals with the decom- 
position of C{X) where X is a compact space on which G acts homogeneously. 
If Gxo is the stabilizer of a base point xq G X, then X can be identified with the 
quotient space G/Gx^- The Haar measure on G gives rise to a G-invariant regular 
measure fi on X and C{X) is endowed with the inner product 

(/,/') ■.= ^Jj{x)J{x)d^,{x). 

The space C{X) can be identified with the space C{G)'^^o of G^p -invariant (for the 
right translation) functions thus C{X) affords a decomposition of the form 



C{X) c^^R^n R 
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for some integers tur, < niR < dR, in the sense of uniform as well as 
convergence. 

A more serious generalization of the above theorem deals with the unitary rep- 
resentations of G. These are the continuous homomorphisms from G to the unitary 
group of a Hilbert space. 

Theorem 3.2. Let tt : G ^ U{H) be a continuous homomorphism from G to 
the unitary group of a Hilbert space H. Then H is a direct Hilbert sum of finite 
dimensional irreducible G-modules. 

Proof. The idea is to construct in H a G-subspace of finite dimension and then to 
iterate with the orthogonal complement of this subspace. To that end, for a fixed 
V & H, one chooses / G C{G) such that Jq f{g){7r{g)v)dg / 0. From Peter 
Weyl theorem, / can be assumed to be a finite linear combination of the Rij. 
In other words, there exists a finite dimensional unitary representation {p, V) and 
ei, 62 G V such that f{g) = {p{g~^)ei,e2)v- The operator T : V ^ H defined 
by 

T(x) = / {pig-')x,e2)v{7^i9)v)dg 
Jg 

commutes with the actions of G and is non zero. Thus its image is a non zero 
G-subspace of finite dimension of H. 

□ 

3.3. Examples. 

3.3.1. The unit sphere S"""^: it is the basic example. The orthogonal group G = 
0(M") acts homogeneously on S"'~^. The stabilizer Gxq of xq can be identified 
with 0{xq ) ~ 0(M"^^). Here = a; is the Lebesgue measure on S'^~^. We set 
coji '■= 00 {S'""^). The irreducible decomposition of C{S'^~^) is as follows: 

where is isomorphic to the space Harm^ of harmonic polynomials: 
Harm^ := {P e C[Xi, . . . , : AP = 0, A = ^ ^} 

i=l * 

The space Harm^ is a 0(M") -module because the Laplace operator A commutes 
with the action of the orthogonal group and it is moreover irreducible. Its dimen- 
sion equals h]^ := ("+^~^) - The embedding of Harm^ into is 
the obvious one, to the corresponding polynomial function in the n coordinates. 

3.3.2. The action of stabilizers of many points: for our purposes we are inter- 
ested in the decomposition of some spaces C{X), X homogeneous for G, for the 
action of a subgroup H of G, typically H = Gxi,...,xs the stabilizer of s points. In 
order to describe it, it is enough to study the decomposition of the G-irreducible 
submodules of C{X) under the action of H; thus we have to decompose only fi- 
nite dimensional spaces. However, because the same irreducible representation 
of H may occur in infinitely many of the G-isotypic subspaces, it happens that 
the iJ-isotypic subspaces are not of finite dimension. A typical example is given 
by X = G = 0(M") and = Ge ~ 0(M"-i). It is a classical result 
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that for the restricted action to H the decomposition of Harm^ into i7-irreducible 
subspaces is given by: 



n-l 



(8) Harm^ ~ ^ Harm 

Hence, each of the in (I3.3.1I ) decomposes likewise: 

Trn TTTi I Trn i i Trn 

-"fc — ^0,k -L J^l,k -L • • • -L J^k,k 

where i/"^ ~ Harm"^^. We have the following picture, where the iJ-isotypic 
components appear to be the rows of the second decomposition. 

=G 1- ± ... ± HI ± . . . 

=H H^fl -L -^0,1 -L ••• -L H^^k 

± HI, ± ... ± HI, 



4. Harmonic analysis of compact spaces 

We take notations for the rest of the lecture notes. X is a compact space (possi- 
bly finite) on which a compact group (possibly finite) G acts continuously. More- 
over, X is endowed with a G-invariant Borel regular measure n for which ijl{X) 
is finite. If X itself is finite, the topology is the discrete topology and the measure 
is the counting measure. In the previous sections we have discussed the decom- 
position of the permutation representation C{X). In order to lighten the notations, 
we assume that G has a countable number of finite dimensional irreducible rep- 
resentations (it is the case if G is a group of matrices over the reals since then 
LP'{G) is a separable Hilbert space), and we let TZ = {Rt, k > 0}, where Rq is 
the trivial representation. We let dk := dim(i?fc). From Theorem 13.21 we have a 
decomposition 

(9) C{X) C L^{X) = ek>0,l<i<m,Hk,^ 

where Hk^i C C{X), Hk^i ~ Rk,0 < nik < +oo (the case = means that Rk 
does not occur, the case nik = +oo may occur if G is not transitive on X). The 
isotypic subspaces are pairwise orthogonal and denoted Ik '■ 

Ik = ®i=lHk,i 

We take the subspaces Hk^i to be also pairwise orthogonal. For all k, i, we choose 
an orthonormal basis Ck^i^i, ■ ■ ■ , ek,i,dj, of Hk^i such that in this basis the action of 
g £ Gis expressed by the unitary matrix Rk{g)- The set {e^ j ^} is an orthonormal 
basis in the Hilbert sense. 

4.1. Commuting endomorphisms and zonal matrices. In this subsection we 
want to give more information on the algebra EndG(C(X)) of commuting con- 
tinuous endomorphisms of C{X). We introduce, for K G C(X^), the operators 
Tk, called Hilbert-Schmidt operators: 

TkU){x) = I^K{x,y)f{y)di,{y). 

It is easy to verify that Tk G EndGr(C(X)) if K is G-invariant, i.e. if K{gx, gy) = 
K{x,y) for all g £ G, {x,y) G X"^. A continuous function K{x,y) with this 
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property is also called a zonal function. It is also easy, but worth to notice that 
Tk o Tx' = Tx*K' where K * K' is the convolution of K and K': 

iK*K'){x,y) := [ K{x, z)K' {z,y)dfi{z). 
Jx 

Let 

/C := {K G C{X^) : K{gx, gy) = K{x, y) for all g^G, (x, y) G X^]. 

The triple (/C, +, *) is a C-algebra (indeed a C*-algebra, with K*{x, y) := K{y, x)). 
Thus we have an embedding /C — EndG(C(X)). 

Assume V C C(X) is a finite dimensional G-subspace such that V = Wi _L 
• • • _L Wm with Wi R € TZ. By the same proof as the one of Theorem 12.41 
EndG(l^) — £mxm_ Morc precisely, we have seen that, if j : Wi —?- Wj are 
G-isomorphisms, such that Ukj o Uj^i = Uk,i, then an element (p G EndG(y) 
is associated to a matrix A = (aij) G ([^rnxm ^^^.j^ ^j^^^^ f^j. ^jj j ^ WiXh 

pwAf) = 

m 

For all 1 < i < m, let (ej^i, . . . , Cj^d), d = dim(i?), be an orthonormal basis of Wi 
such that in this basis the action of 5 G G is expressed by the unitary matrix R{g). 
We define 

d 

Ei,j{x,y) ■.= ^ei^s{x)ej^s{y)- 

s=l 

Then we have: 

Lemma 4.1. The above defined functions Eij satisfy: 

(1) Eij is zonal: Eij{gx,gy) = Eij{x,y). 

(2) Let Tij := Te,^^. Then Tj^Wi) = Wj and Tj^Wk) = 0/or k ^ i. 

(3) Tij o Tj^k = Ti^k- 

Proof. (1) From the construction, we have 

d 

ei,s{gx) = y^Rs,t{g)ei,t{x) 
t=i 

thus 

d 

Ei,j{gx,gy) = ^ei^s{gx)ej^s{gy) 

s=l 
d d 

s=l k,l=l 
d d 

= XI {j^^s,k{9)Rs,i{9))^iA^yjAy) 

k,l=l s=l 
d 

= ^ei^k{x)ej^k{y) = Eij{x,y) 

k 
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where the second last equality holds because R{g) is a unitary matrix. 

(2) We compute Tj^i{ek,t)- 

Tj,i{ek,t){x) = J^{Ylej,s{x)ei,s{y))ek,t{y)dn{y) 

d 

d 
s=l 

(3) Similarly one computes that 

□ 

The Eij{x, y) put together form a matrix E = E{x, y), that we call the zonal 
matrix associated to the G-subspace V: 

(10) E{x,y):={E,,{x,y))^^^^^. 

At this stage is is natural to discuss the dependence of this matrix on the various 
ingredients needed for its definition. 

Lemma 4.2. We have 

(1) E{x^ y) is unchanged if another orthonormal basis ofWi is chosen (i.e. if 
another unitary representative of the irreducible representation R is cho- 
sen). 

(2) E{x, y) is changed to AE(x, y)A* for some matrix A G G1(C™) if another 
decomposition (not necessarily with orthogonal spaces) V = W[ © • • • ® 

is chosen. 

Proof. (1) Let (e^ i, • • • , ^) be another orthonormal basis of Wj and let Ui 
be unitary dx d matrices such that 

= (ej,i,---,ej,d)C/i- 

Since we want the representation R to be realized by the same matrices in 
the basis (e' . . . , ^) when i varies, we have Ui = Uj = U. Then, with 
obvious notations, 

E'i,j{x, y) =(e-^i(x), . . . , e-^^(a;))(e- . . . , 6-^^(2/))* 

=(ej,i(x), . . . ,ei^d{x))UU*{ei^i{y), Ci^diy))* 
=(ei,i(x), . . . , ei^d{x)){ei,i{y), • • • , ei,d{y)T 
=Eij{x,y). 

(2) IfV = Wi ± ••• J.Wm = Wil. ••• ± W-;;, with basis (ei,i,...,ei,d) 
of Wi and (e^ 1, • • • , ^) of PF/ in which the action of G is by the same 
matrices R{g), let (f) G End(y) be defined by (t>{ei^s) = s- Clearly 
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(p commutes with the action of G; if Uj^i is defined by Uj^i{ei^s) = ej^s 
then we have seen that, for some matrix A = (aij), ^ = 4>{ei^s) = 
SjLi '^j,i^j,s- Moreover A is invertible. It is unitary if the spaces Wl are 
pairwise orthogonal. With the notations E{x) := (ej^s(x)), we have 

E{x,y) = E{x)E{y)* and E'{x) = A*E{x) 

thus 

E'{x,y) = A'E{x,y)A. 

□ 

Going back to (/> G Endcl^), from Lemma |4~T] we can take Uj^i = Tj j and we 
have the expression 

m 

We take the following notation: the space of linear combinations of elements of 
the form f{x)g{y) for (/, g) G V'^ is denoted V^'^\ We have proved the following: 

Proposition 4.3. Let 

Kv := {K eV^^^ ■.K{gx,gy) = K{x,y)forallgeG, {x,y)eX^]. 

The following are isomorphisms ofC-algebras: 

ICv ^EndG(V^) C"^™ ^EndG(l^) 

K ^Tk a ^T^a,e)- 

Moreover, EndG(C(X)) is commutative iff fC is commutative iffnik = 0, 1/or all 
k>0. 

Proof. The isomorphisms are clear from previous discussion. For the last assertion, 
it is enough to point out that 

EndG(C(X)) = J]EndG(Xfc). 

A;>0 

□ 

Remark 4.4. Proposition \4.3\ shows in particular that ICv and FindclV) have 
the same dimension. It is sometimes easy to calculate the dimension of fCy; for 
example if X is a finite set and V = C{X), then dim(/Cy) is exactly equal to 
the number of orbits of G acting on X^. On the other hand, in this case, the 
dimension ofEjudciV) is the sum of the squares of the multiplicities in C{X) of 
the irreducible representations of G. For the binary Hamming space treated in 
\2.8.1\ the orbits of G acting on X"^ are in one to one correspondance with the 
values taken by the Hamming distance, i.e. there are (n + 1) such orbits. Thus, 
once we have obtained the decomposition C{Hn) = Pq -L ■ ■ ■ 1. Pn, because this 
decomposition involves allready (n + 1) subspaces, we can conclude readily that 
these subspaces are irreducible. This reasoning applies also to the Johnson space 
\2.8.2\ and to the more general q-Hamming space \573\ A variant of this method is as 
follows: if we suspect V C C{X) to be irreducible, then it is enough to prove that 
ICv has dimension 1. See in \5.3.1\ for an illustration. 
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4.2. Examples: G-symmetric spaces. 

Definition 4.5. We say that X is G-symmetric if for all (x, y) G X'^, there exists 
g £ G such that gx = y and gy = x. In other words, {x, y) and {y, x) belong to 
the same orbit of G acting on X^. 

A first consequence of Proposition 14.31 is tliat G-symmetric spaces liave multi- 
plicity free decompositions. 

Proposition 4.6. If X is G-symmetric then rrik = 0, 1/or all k > and Ek{x, y) 
is real symmetric. 

Proof. For all K £ IC, K{x, y) = K{y, x). Thus K, is commutative: indeed, 



4.2.1. 2-point homogeneous spaces: these spaces are prominent examples of G- 
symmetric spaces. 

Definition 4.7. A metric spaces (X, d) is said to be 2-point homogeneous for the 
action of G if G is transitive on X, leaves the distance d invariant, and if, for 



Examples of such spaces of interest in coding theory are numerous: the Ham- 
ming and Johnson spaces, endowed with the Hamming distance, for the action of 
respectively T ys Sn and Sn', the unit sphere S*"^^ for the angular distance 6{x, y) 
and the action of the orthogonal group. It is a classical result that, apart from 
the projective spaces P"(A') for K = M, C, M, and p2(o), are the only real 
compact 2-point homogeneous spaces. 

There are more examples of finite 2-point homogeneous spaces, we can mention 
among them the g- Johnson spaces. The g- Johnson space Jn{q) is the set of linear 
subspaces of of fixed dimension w, with the action of the group Gl(Fg ) and the 
distance d{x, y) = dim(x + y) — dim(3; n y). We come back to this space in the 
next section. 

There are other symmetric spaces occurring in coding theory: 

4.2.2. The Grassmann spaces: X = Qm,n{K), K = R,C, i.e. the set of m- 
dimensional linear subspaces of K", with the homogeneous action of G = 0(M") 
(respectively [/(C")). This space is G-symmetric but not 2-point homogeneous (if 
m > 2). The orbits of G acting on pairs {p, q) G X"^ are characterized by their 
principal angles [21]. The principal angles of {p, q) are m angles (^i, . . . , 6m) G 
[0,7r/2]'" constructed as follows: one iteratively constructs an orthonormal basis 
(ei, . . . , em) of p and an orthonormal basis (/i, . . . , fm) of q such that, for 1 < 




(K * K'){y, x) = {K* K'){x, y). 



Moreover Ek{x, y) = Ek{x, y) 



Ek{y,x). 



□ 
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i < m, 



cos 6i 



max{|(e,/)| 



eGp, f £q, 
(e,e) = (/,/) = !, 

(e,e,) = (/,/,) = fori <j <i-l} 



= \iei,fi)\ 
The we have (see |[211 ): 



there exists g £ G such that ((7p, (75) = {p', q') 
ie.ip, q),..., aUP, q)) = (Oiip', q'),---, Om{p', q'))- 



4.2.3. The ordered Hamming space: X = (Fg)" (for the sake of simphcity 
we restrict here to the binary case). Let x = (xi, . . . ,x„) G X with Xi G Fj. 
For y G F2, the ordered weight of y, denoted Wr{y), is the right most non zero 
coordinate of y. The ordered weight of x G X is Wr{x) := X^iLi ^^(^^i) and the 
ordered distance of two elements (x, y) G X"^ is dr{x, y) = Wr{x — y). Moreover 
we define the shape of (x, y): 



Another expression of Wr{x) is Wr{x) = iej. 

If B is the group of upper triangular matrices in G1(F2), and i^aff the group of 
affine transformations of F2 combining the translations by elements of F2 with B, 
the group G := xi Sn acts transitively on X. Since B acting on Fg leaves Wr 
invariant, it is clear that the action of G on X leaves the shape shape(x, y) invari- 
ant. More precisely, the orbits of B on Fg are the sets {y G Fg : Wr{x) = i} and, 
consequently, the orbits of G acting on X"^ are characterized by the so-called shape 
of (x, y). Since obviously shape(x, y) = shape(y, x) it is a symmetric space. This 
space shares many common features with the Grassmann spaces, especially from 
the point of view of the linear programming method (see Q, jOl, IIBTI ). 

4.2.4. The space X = T under the action of G = F x F: the action of G is by 
(7, 7')x = 7x7'"^. Then two pairs (x, y) and (x', y') are in the same orbit under 
the action of G iff xy~^ and x'y'^^ are in the same conjugacy class of F. Obviously 
(x,y) and {y~^,x~^) are in the same G-orbit. We are not quite in the case of a 
G-symmetric space however the proof of the commutativity of /C of Proposition 
14. 61 remains valid because the variable change x — )• x^^ leaves the Haar measure 
invariant. 

4.3. Positive definite functions and Bochner theorem. 

Definition 4.8. A positive definite continuous function on X is a function F G 
C(X^) such that F{x, y) = F(y, x) and one of the following equivalent properties 




hold: 



(1) For all n, for all (xi , . . . , x„) G X'"-, for all {ai , . . . , q„) G C^, 



n 
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(2) For all a & C{X), 



a{x)F{x,y)a{y)dn{x,y) > 0. 

This property will be denoted F ^ 0. 

The first property means in otlier words tliat, for all choice of a finite set of 
points {xi, . . . ,Xn) S X^, the matrix {F(xi,Xj))i<,ij<n is hermitian positive 
semidefinite. The equivalence of the two properties results from compactness of 
X. Note that, if X is finite, F is positive definite iff the matrix indexed by X, with 
coefficients F{x, y), is positive semidefinite. 

We want to characterize those functions which are G-invariant. This charac- 
terization is provided by Bochner in 111] in the case when the space X is G- 
homogeneous. It is clear that the construction of previous subsection provides 
positive definite functions. Indeed, 

Lemma 4.9. if A ^ 0, then {A, E) is a G-invariant positive definite function. 
Proof. Let a{x) G C{X). We compute 



a 

X2 



„ m 

{x){A,E)a{y)dfM{x,y) = / Aija{x)Eij{x,y)a{y)dfj,{x,y) 

^^^^ 

m „ 

/ a{x)Eij{x,y)a{y)dfi{x,y) 



m d 



y^^i-i / ci:i^)(ii,s{x)ej,s{y)a{y)d^i{x,y) 

U=ls=l 



m d 

i,j=l s=l 
d m 

s=l i,j=l 

where aj^s := fx (^ix)^i,six)dfi{x). □ 

Remark 4.10. The following properties are equivalent, for a m x m matrix func- 
tion E{x, y): 

(1) For all At 0, {A,E{x,y)) h 

(2) For all (xi, . . . G X", (oi, G C", aiE{xi,Xj)a] h 0. 

r/je proof is left to the reader as an exercise (hint: use the fact that the cone of 
positive semidefinite matrices is self dual). 

To start with, we extend the notations of the previous subsection. We define 
matrices E^ = Ek{x, y) associated to each isotypic component Z^, of size nik x nik 
(thus possibly of infinite size) with coefficients Ek^ij{x, y) defined by: 

dk 

Ek,i,jix,y) ■■= y^^ek,i,six)ek,j,siy)- 

s=l 
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If Fk = {fk,i,j)i<i,j<mk is hermitian, and if J^ij < +00, the sum 

{Fk,Ek) := fk,i,jEk,i,j 

is L^-convergent since the elements ek,i,s{x)ei.j^t{y) form a complete system of 
orthonormal elements of C{X'^). We say is positive semidefinite {F^ ^ 0) if 
^ifk,i,j^j > for all (Ai)i<i<mfc such that ^ |Aip < +00. Then, with the 
same proof as the one of Lemma 1491 the function {F^, E^) is positive definite if 
Fk h 0. The following theorem provides a converse statement (see UJJ). 

Theorem 4.11. F G C{X'^) is a G-invariant positive definite function if and only 
if 

(11) F{x,y) = Y,{Fk,Ek{x,y)) 

k>0 

where, for all k >0, 

= 1 — TTTfT / F{x,y)Ek{x,y)di2{x,y) h 0, 

and the sum converges to F for the topology. If moreover G acts homogeneously 
on X, the sum (1111) itself converges uniformly. 



Proof. The elements j s(3;)e/ j form a complete system of orthonormal ele- 
ments of C{X^). Hence F has a decomposition 

F{x,y) = J2 fk,i,s,i,j,tek,iA^)ei,j,t{y) 

where the convergence of the sum is L^. The condition F{gx,gy) = F{x,y) 
translates to: 

fk ,i,u,l,j,v — fk,i,s,ij,tRk,u,sig)Ri,v,tig)- 

s.t 

Integrating on G G and applying the orthogonality relations of Theorem 12.5 
shows that fk,i,u,i,j,v = if / Z or u / Moreover it shows that fk,i,u,k,j,u does 
not depend on u. The resulting expression of F reads: 

F{x,y) = ^ (^^fk,i,jEkAjix,y) 

k>0 i,j 

and 

dkfk,i,j = ^(^2) J^^ F{x,y)Ek,ij{x,y)dfi{x,y), 

which is the wanted expression, with Fk := {fk,i,j)i<i,j<m^.- 

Now we show that Fk >i 0. Let, for k,s fixed, a{x) = J2i'^i^k,i,s{x), with 
Y2i < +00. By density, property (2) of Definition 14. 8 1 holds for a G L^(X). 
We compute like in the proof of Lemma |49] 

"ifc 

a{x)F(x,y)a{y)d^{x,y) = aifk,i,jaj 

thus Fk h 0. 

In the case of X being G-homogeneous, the uniform convergence of the sum in 
(fTTI ) is proved in till- 
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□ 



In order to reduce linear programs involving G-invariant positive definite func- 
tions to finite dimensional semidefinite programs, we need to be able to approxi- 
mate such functions uniformly with finite sums of the type (ITT]) , in other words by 
functions built form finite dimensional subspaces of C{X). A necessary condition 
is thus that all continuous functions on X are uniformly approximated by elements 
of some sequence of finite dimensional subspaces of C{X). Such subspaces are 
usually provided by the polynomial functions of bounded degree, when it makes 
sense. More generally, let us assume that there exists a sequence {Vd)d>o of finite 
dimensional G-subspaces of C{X) such that C Vd-^-i, and Ud>oVd is dense in 
C{X) for the topology of uniform convergence. For example, Peter- Weyl theo- 
rem provides such subspaces when X is F-homogeneous, for a compact group F 
containing G. Then we have the following result: 

Theorem 4.12. Under the above assumptions, if moreover X is homogeneous un- 
der a larger compact group T, and if the irreducible subspaces H^^i are chosen so 
that Hk^i C Vdfor all 1 < i < rrid^k where md^k the multiplicity of Rk in Vd, 
then a G -invariant positive definite function F G C(X^) is the uniform limit of a 
sequence of positive definite functions Fd ^Vd®Vd thus of the form 



where Fd^k ^ matrix of size rtid^k (and thus the sum has a finite number ofnon 
zero terms). 

Proof. We proceed like in the proof of Peter Weyl theorem. Compact self-adjoint 
Hilbert-Schmidt operators on C(X^) are of the form 



We start to construct K such Tk{F) >z and \\Tk(F) — F\\^ is arbitrary small. 
The first condition is fulfilled if K can be expressed in the form K{{x, y), (z, t)) = 
Kq{x, z)KQ{y, t) where Kq{x, z) = Kq{z, x). We take (pQ a continuous function 
on F; if (pQ denotes the left and right average of 4>q over Fq (where X = F/Fq), 
we take Kq{x, y) = 4>'q{'^~^5) for any 7 G x, 5 € y). Then with a suitable choice 
of 00, \\Tk{F) — F\\ao < e (thanks to uniform continuity of F, it is enough that 
(f)Q has support contained in some prescribed open neighborhood of 1, takes values 
between and 1, satisfies (/>o(7) = </'o(7^^) and /p 0o = iTol)- Moreover, Kq is 
F-invariant. 

We can find d > and Lq{x, y) & Vd® Vd such that Lq{x, z) = Lo{z, x) and 
II-Lq — i^olloo is arbitrary small. Replacing Lq by its average on G will not change 
these three properties of Lq. Then, if L{{x, y), {z, t)) := Lo{x, z)Lo{y, t), Tl{F) 
comes arbitrary close to Tk{F) for || ||oo and Tl{F) G Vd. Now, Tl{F) >z 0, 
is invariant under G and belongs to the finite dimensional space Vd (8> Vd thus it has 
the announced form from Theorem 14.1 II 



(12) 




k>0 




□ 



Now the main deal is to compute explicitly the matrices Ek{x,y) for a given 
space X. The next section gives explicit examples of such computation. 
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5. Explicit computations of the matrices Ek{x,y) 

We keep the same notations as in previous section. Since the matrices Ek{x, y) 
are G-invariant, their coefficients are functions of the orbits of G acting on X"^. 
So the first task is to describe these orbits. Let us assume that these orbits are 
parametrized by some variables u = (uj). Then we seek for exphcit expressions 
of the form 

Ek{x,y) = Yk{u(x,y)). 

The measure /U induces a measure on the variables that describe these orbits, for 
which the coefficients of are pairwise orthogonal. This property of orthogonal- 
ity turns to be very useful, if not enough, to calculate the matrices E^. 

The easiest case is when the space X is 2-point homogeneous for the action of 
G, because in this case the orbits of pairs are parametrized by a single variable 
t := d{x, y). Moreover we have already seen that in this case, the decomposition 
of C{X) is multiplicity free so the matrices Ek{x, y) have a single coefficient. 

5.1. 2-point homogeneous spaces. We summarize the results we have obtained 
so far: 

C{X) = ®k>oHk 

where are pairwise orthogonal G-irreducible subspaces; to each is associ- 
ated a continuous function Pk{t) such that Ek{x, y) = Pk{d{x, y)) and 

F^O^ F = Y, fkPk{d{x, y)) with fk > 0. 

fc>0 

Pk{t) is called the zonal function associated to Hk- Since the subspaces Hk are 
pairwise orthogonal, the functions Pk{t) are pairwise orthogonal for the induced 
measure. This property of orthogonality is in general enough to determine them in 
a unique way. We can also notice here that Pfc(O) = dk- This value is obtained 
with the integration on X of the formula i-*A:(0) = ^1^=1 ^k,i.s{x)ek.i,s{x)- 

5.2. X = {1, . . . ,g} under the action of Sg. This is a very easy case, which 
will play a role in the study of the g-Hamming space. Since the constant function 
1 is 5q-invariant, we have the Sg decomposition C{X) = C 1 _L L. Obviously, 
the action of Sg on X^ has two orbits: the set of pairs {i, i), and the set of pairs 
(i, j) for i 7^ j. Thus, from Proposition 14.31 and Remark l44l L is irreducible. We 
let zq := 1 and choose an orthonormal basis (zi, . . . , Zg-i) of L. We want to 
compute the zonal function El associated to L. We have by definition El{x, y) = 
Y11z\ Zi{x)zi{y) and El takes only two different values: one for x = y and one 
for X ^ y. We have El{0, 0) = dim(L) = q — 1 and we can compute El{0, 1) 
easily using the fact that £'l(0, y) is orthogonal to zq thus J2l=i El{0, u) = = 
El{0, 0) + {q- 1)El{0, 1). Thus ^^(0, 1) = -1. 

5.3. The g-Hamming space. In the binary case we have already calculated the 
functions Pk{t) in l2.8.11 Indeed, the irreducible subspaces Pk afford the orthonor- 
mal basis {xz,wt{z) = k}. So, 

Ek{x,y)= Xz{x)Xz{y)= E i-ir'^'^^'^ =Kk{dH{x,y)) 

wt{z)=k wt{z)=k 

from ([Hi. Now we treat the more general g-Hamming space. This is the space 
Hn,q = where F is a finite set with q elements denoted F = {ao, ai, . . . , ag-i}. 
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The semidirect product G = S'g x 5„ acts on Hn,q and leaves the Hamming dis- 
tance invariant. Here the permutation group Sq acts on F by raj = a^-i^i^ while the 
permutation group 5„ acts on Hn,q by a{xi, . . . , x„) = {x^-i(^i^, . . . , x„-i(^n))- 
Moreover G acts on Hn,q 2-point homogeneously. The action of Sq on C{F) is 
studied in 15.21 and we take the same notations. We define = (^i, . . . , (/>„) G 
C{Hn,q) where e {zq, zi,..., Zq-i} by: = nr=i 4>i{xi)- These elements 
(j) form an orthonormal system: it is easy to see that 

n 
i=l 

We define the weight of (p by: wt{cl)) := \{l < i < n : (/){ zo}\. For < A; < n, 
let Pk be the subspace generated by the set of (j) with wt{4>) = k. The dimension 
of Pk is the number of such (f>, which is equal to {q — l)'^(^) and we have the 
decomposition 

(13) CiHn) =Po±Pl±--- ±Pn. 

An element t £ Sq act trivially on zq and sends Zj for i 7^ to a linear combination 
of zi, . . . , Zg-i. Thus for all g £ G, g(p is a. linear combination of ^'s with the 
same weight as (p and G stabilizes P^. The action of G on pairs of elements of 
Hn,q has exactly {n + 1) orbits corresponding to the {n + 1) values 0, 1, . . . , n 
that the Hamming distance takes thus we can conclude that Pk is irreducible from 
Proposition 14.31 Now we compute the zonal function Ek{x, y) attached to Pk- By 
definition we have 

Ek{x,y)= ^ (l}{x)4>{y) 

(j),wt{4')=k 

and we want to calculate P^ such that Pk{t) = Ek{x,y) for any {x,y) with 
d{x, y) = t. We set a; = (ai, . . . , ai, ao, • • • , ao) where t coordinates of x are 
equal to ai and y = (ao, • • • , ao). For all cp, we let i := |{j : 1 < j < n : Xj = 
ai and (pj zq} and reorder the set of (p € Pk according to i. 

Pk{t) = r j K _ J 5^ n ^i«(«i)^jV.(«o) n ^jV,("o)^jJao) 
j=o ^ ^ ^ ii,- -,ifc^o?i=i M=i+i 

i-n V / \ / s=l s=l ' 

I -El (ai , ao ) ' (ao , ao ) 



=0 

k 

E 



t\/n-t 
i J \k — i 

t\ f n — t 



k 



^ \i J \ k — i 

i=0 



-inq-1)' 



with the notations and results of 15.21 Pk{t) is equal to the Krawtchouck polyno- 
mial K^''^{t) of parameters q and n which satisfies the following characteristic 
properties: 

(1) deg(i^,"'^) = k 

(2) i^no) = ('?-i)'a) 
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(3) Orthogonality relations: for all < /c < / < n 

The orthogonality relations aie direct consequences of the orthogonality of the 
subspaces Pfc. 

5.3.1. The Johnson space J^: with the notations of subsection 12.8.21 we have 
shown the decomposition 

C{J:^) ^ ± H^-i ± ■ ■ ■ ± Ho 

but not yet the irreducibility of Hi. So far their might by several Pij, j = 1, . . . 
associated to H^. The zonal functions express as functions of t := |x n y| the 
number of common ones in x and y. The orthogonality relation is easy to compute: 

n 

Yl fi\xr\y\)f'i\xny\) = ^card{x : \x H y\ = i} f 

x€X i=0 



± fin, -i)fW^). 



i=0 

By induction on k one proves that P^j has degree at most k in t. The conditions: 

(1) degiQk) = k 

(2) Qfc(O) = 1 

(3) for all < A; < Z < n 



g(T)("r)*««'«=» 



determine a unique sequence {Qq, Qi, . . . , Qw)- Thus there is only one Pj^ j for 
each k and it is equal to hkQi.{w — t). The polynomials Qk defined above belong 
to the family of Hahn polynomials. 

5.3.2. The sphere 5"^^ : the distance on the sphere is the angular distance 6{x, y). 
It appears more convenient to express the functions in the variable t = x ■ y = 
cos 6{x, y). A standard calculation shows that 

f{x ■ y)df,{y) = Cn / f{t){l - t^)^dt 

for some irrelevant constant c^. The conditions: 
. deg(Pfc") = k 

• ^'fcni) = 1 

• For all ki^l, j\ P^{t)Pp{t){l - t^)'^dt = 

define a unique sequence of polynomials by standard arguments (i.e. obtained by 
Gram Schmidt orthogonalization of the basis (1, . . . , t'^, . . . )), it is the sequence 
of so-called Gegenbauer polynomials with parameter n/2 — 1 B3l . The decompo- 
sition [33TT]of C(S'"~^) shows that, to each /c > the function Pk{x ■ y) associated 
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to HJ^ ~ Harm^ is polynomial in x ■ y and satisfies the above conditions except 
the normalization of i-fc(l) thus we have Pkit) = h^P^{t). 

5.3.3. Other 2-point homogeneous spaces: as it is shown in the above exam- 
ples, a sequence of orthogonal polynomials in one variable is associated to each 
such space. In the case of the projective spaces, it is a sequence of Jacobi polyno- 
mials. We refer to |[24ll . |[28l . BSl for their determination in many cases and for the 
applications to coding theory. 

5.4. Other symmetric spaces. Now we turn to other cases of interest in coding 
theory, where the space X is symmetric but not necessarily 2-point homogeneous. 
Since the decomposition of C{X) is multiplicity free, the matrices Ek{x,y) still 
have a single coefficient which is a member of a sequence of orthogonal polynomi- 
als, but this time multivariate. The first case ever studied (at least to my knowledge) 
is the case of the non binary Johnson spaces 1441 . its associated functions are two 
variables polynomials, a mixture of Hahn and Eberlein polynomials. We briefly 
discuss a few of these cases. 

5.4. 1. The Grassmann spaces: Q the orbits of X'^ are parametrized by the prin- 
cipal angles (^i, . . . , 9m) (14.2.21 ). The appropriate variables are the yi := cos^ 6i. 
The decomposition of C{Qm,n) under 0(M") (respectively [/(C")) together with 
the computation of the corresponding sequence of orthogonal polynomials was 
performed in ll23l . We focus here on the real case. We recall that the irreducible 
representations of 0(M") are (up to a power of the determinant) naturally indexed 
by partitions k = . . . , where ki > • • • > k„ > (we may omit the 
last parts if they are equal to 0). Following |[22]| . let them be denoted by F^. For 

(k) 

example, Vn = C 1, and Vn = Harrrifc. 

The length ^(k) of a partition k is the number of its non zero parts, and its degree 
deg(K) also denoted by |«;| equals Yli=i ^i- 

Then, the decomposition of C{Qm,n) is as follows: 

where k runs over the partitions of length at most m and 2k stands for parti- 
tions with even parts. We denote by P^iyi, - ■ ■ ,ym) the zonal function associ- 
ated to y^''. It turns out that the are symmetric polynomials in the m variables 
yi , . . . , ym, of degree | k| , with rational coefficients once they are normalized by the 
condition Pk{1, • • • , 1) = 1- Moreover, the set {PK)\K\<k is a basis of the space 
of symmetric polynomials in the variables yi, ... ,ym of degree at most equal to k, 
which is orthogonal for the induced inner product calculated in Il23l . 

m m 

i,j=l i=l 
i<j 

(One recognizes a special case of the orthogonal measure associated to generalized 
Jacobi polynomials ( |[25]| ). 
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5.4.2. The ordered Hamming space: it follows from the discussion in l4.2.3l that 
the variables of the zonal functions are the (cq, ei, . . . , e,.). Elaborating on the 
computation explained above for the Johnson space, one can see that in the case 
of finite spaces, the weights of the induced measure are given by the number of 
elements of the orbits of X under the action of Stab(e) for any e S X. Taking 
e = 0^", thus Stab(e) = xi 5„, and the orbit of x is the set of elements with the 
same shape (/o, . . . , fr) as x. The number of such elements is (j^^j )2^i(*~^)'^\ 
These are the weights associated to the multivariate Krawtchouk polynomials. 

5.4.3. The space X = T under the action of G = F x F: we need an ex- 
plicit parametrization of the conjugacy classes of F, which is afforded by very few 
groups. Famous examples (if not the only ones) are provided by the permutation 
groups and the unitary groups. In the first case the parametrization is by the de- 
composition in disjoint cycles and in the second case it is by the eigenvalues. The 
decomposition ofC{X) is given by Peter Weyl theorem 

C(F) = ^R^R* 

Ren 

and the associated functions Pr{x, y) are the characters: 

PR{x,y) = XR{xy'^)- 

In both cases (5„ and [/(C")) the irreducible representations are indexed by par- 
titions A and there are explicit expressions for P\. In the case of the unitary 
group P\{xy^^) are the so-called Schur polynomials evaluated at the eigenvalues 
of xy~^. 

5.5. Three cases with non trivial multiplicities. So far the computation of the 
matrices Ef^{x, y) in cases of non trivial multiplicities has been worked out in very 
few cases. We shall discuss three very similar cases, namely the unit sphere of 
the Euclidean sphere ([4|), the Hamming space (1.461), and the projective geom- 
etry over Fg (17]), where the group considered is the stabilizer of one point. In 
the case of the Hamming space, this computation amounts to the computation of 
the Terwilliger algebra of the association scheme and was performed initially by 
A. Schrijver in BOl . who treated also the non binary Hamming space |[20l . The 
framework of group representations was used in ll46l to obtain the semidefinite 
matrices of 1*401 in terms of orthogonal polynomials. We present here the uniform 
treatment of the Hamming space and of the projective geometry in the spirit of ifTTl 
adopted in ||7l. We also generalize to the case of the stabilizer of many points in 
the spherical case and enlighten the connection with the positive definite functions 
calculated in fM\. 

5.5.1. The unit sphere 5""^ with G := Stab(e, 0(M'')). We continue the dis- 
cussion initiated in l3.3.2l and we follow [4|. Let E^{x, y) be the zonal matrix as- 
sociated to the isotypic subspace related to Harm^~^ and to its decomposition 
described in 13.3.21 

1-k = Hk,k -L Hk,k+1 -L • • • 

We index with i,j > so that E'^^-{x^y) is related to the spaces HJ}^_^-, 
H^k+j- "^^^ orbits of G on pairs of points {x, y) G X'^ are characterized by the 
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values of the three inner products u := e - x,v := e-y and t := x-y. Thus (li, v, t) 
are the variables of the zonal matrices and we let: 

El{x,y) = Y^{u,v,t). 

Theorem 5.1. [^] 

(14) = \kAk,3P?^^\u)P^+^\v)Ql-\u,v,t), 

where 

and Afc j are some real constants. 

Proof. We need an explicit construction of the spaces H^^^^- We refer to HJ 
Ch. 9.8]. For x £ 5""^ let 

X = ue+ \/l — u'^C, 

where u = x-e and C belongs to the unit sphere S'"-^ of (Me)-^. With / G HJ^'^ C 
associate v3(/) G C(5'^-i) defined by: 

v.(/)(x) = (i-^Y/V(C). 

Moreover, we recall that HJ^ is a subspace of the space Pol<fc(S"^~^) of polynomial 
functions in the coordinates of degree at most k. Note that the multiplication by 
(1 — n^)^/^ forces (p{f) to be a polynomial function in the coordinates of x. Clearly 
(p commutes with the action of G. Hence ip{H'^~^) is a subspace of Pol<fc(5"'~^) 
which is isomorphic to Harm^"^. It is clear that these spaces are pairwise or- 
thogonal. More generally, the set {(p{f)P{u) : / G Harm^~^, deg P < i} is a 
subspace of Pol<;j+j(S'"^^) which is isomorphic to i + 1 copies of Harm^^^. By 
induction on k and i there exist polynomials Pi{u) of degree i such that -f^^^^j := 
(p{Hj^^^)Pi{u) is a subspace of HJ^_^_^. This construction proves the decomposition 
dSJ. Moreover, we can exploit the fact that the subspaces Hj^J^ are pairwise or- 
thogonal to prove an orthogonality relation between the polynomials P^. Then this 
orthogonality relation will enable us to identify the polynomials Pi with Gegen- 
bauer polynomials, up to the multiplication by a constant factor. Let us recall that 
the measures on S"^^ and on 5""^^ are related by: 

dcOnix) = (1 - u2)('^-3)/2dud6<J„_i(C). 

Whenever i ^ j we have for all / G HJ^~^ 

= — / ip{f)Pi{uMf)Pjiu)dUn{x) 

= — [ \fiO\Hl - u')''P,{u)P^)du^{x) 

= — [ |/(C)Pd^n-l(C) /' (1 - tx2)'=+("-3)/2p^(^)7^d^, 
J5"-2 J-i 

from which we derive that 

(1 - u2)'^+("-3)/2p.(^)7r(^dn = 0; 
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hence the polynomials Pi{u) are proportional to P" (u) (thus with real coef- 

rn— 1 



ficients..)- We obtain an orthonormal basis of Hj^ from an orthonormal basis 

<n 



(/i, . . . , A) of H'^ 1 by taking ek,i,s = \kMfs)P?'^'^''{u) for a suitable normal- 
izing factor Afc j > 0. With these basis we can compute E^- : 



Ek,ijix,y) = ek,t,s{x)ek,j,s{v) 

s=l 



s=l 

where we have written y = ve + VT^-^C and where the last equality results from 
the analysis of zonal functions of 5"""^. Since 



we have completed the proof. □ 

5.5.2. The unit sphere 5""^ with the action of G := Stab(ei, . . . , e,, 0(M")). 
We assume that (ei, . . . ,6^) is a set of orthonormal vectors. The group G := 
Stab(ei, . . . , Cs, 0(M")) is isomorphic to 0(]R"^'^). The orbit of a pair (x, y) G 
under G is characterized by the data: t := x ■ y, u := {x ■ ei, . . . , x ■ eg), v := 
(y • ei, . . . , y • e^). The decomposition ([8]l applied recursively shows that C(5'"~^) 
decomposes as the sum of G-iiTcducible subspaces where k = {ko, . . . , kg), 
ko < ki < ■ ■ ■ < ks, with the properties: 

Hk C C Polfc, , Hk ^ Harm^-" 

where A;*-*"^ = {kg-r+i, ■ ■ ■ ,ks). Thus, for a given k^, the multiplicity of the iso- 
associated to Harm^~ 



typic component associated to Harm? ^ in Pol<d is the number of elements 



of 

Kd ■■= {{ki, . . . ,ks) : ko < ki < ■ ■ ■ < ks < d}. 
We construct the spaces like in the proof of Theorem l5.lt for x £ let 



X = uiei H h UgCs + \/l - \u\'^C 

where u = (wi,...,n,) and \u\'^ = ELi ■ Let ip : ^ C(5"-i) be 

defined by ip{f){x) = (1 - \u\^)''°/'^ f{C). Then (fiH'^^") = H^s+i where = 
{ko, ko, . . . , ko) and we set, for / = (Zi, . . . , Is), Hko,i '■= ■ ■ • u^^H^s+i. It is 

clear that H^g^i Harm^^'' and that H^q^i C Pold if li + ■■■ + ls < d — k^ thus, 
since 

K'^:={l = {h,...,ls) : k>0, h + --- + ls<d-ko} 
has the same number of elements as K^, 

^ko = ®l&K'^Hko,l- 
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This sum is not orthogonal but we can still use it to calculate E'fc,,, the change 
will be to AEk{x, y)A* for some invertible matrix A. The same calculation as in 
Theorem 15 ■ 1 1 shows that, (up to a change to some AYf^A*): 

Yk,ij{u, V, t) = u^-''vl-^Ql-'{u, V, t) 
with the notations: u^~'' := u'^^^u'^^^ . . . ■u*/"'' and 

g«-(n,^,t) = ((1 - |np)(l - |^|2)fc/2pn-s/ t-{u-v) 

With Bochner Theorem l4.11l we recover the description of the multivariate positive 
definite functions on the sphere given in Il34l . 

5.5.3. The Hamming space and the projective geometry. The set of all Fg- 

linear subspaces of F^, also called the projective geometry, is denoted by V{n, q). 
The linear group Gl(n,Fg) acts on V{n, q). The orbits of this action are the sub- 
sets of subspaces of fixed dimension, i.e. the g- Johnson spaces. If the Hamming 
space Fg is considered together with the action of the symmetric group 5„, the 
orbits of this action are the Johnson spaces. In 1 17] the Johnson space and the q- 
Johnson spaces are treated in a uniform way from the point of view of the linear 
programming method, the latter being viewed as g-analogs of the former. Thus the 
Johnson space corresponds to the value g = 1. In particular the zonal polynomials 
are computed and they turn to be g-Hahn polynomials. Here we want to follow the 
same line for the determination of the zonal matrices E{x, y) in both cases. 

We take the following notations: if (7 is a power of a prime number, we let 
X = V{n, q) and G = Gl(n, Fg), and, li q = 1, we let X be the Hamming space, 
identified with the set of subsets of {1, . . . , n}, and G = Sn the symmetric group 
with its standard action on X. Let 

, I _^ f wt{x) ifq = l 
' ' ' I dim(x) iiq> 1 

For all = 0, . . . , n, the space is defined by 

Xw = {x G X : |x| = w}. 

These subsets of X are exactly the orbits of G. The distance on X is given in every 
case by the formula 

(15) d{x,y) = \x\ + \y\ — 2\x r\y\. 

The restriction of the distance d to Xyj equals d{x, y) = 2{w — \x (^y\) and it is 
a well known fact that G acts 2-points homogeneously on X^- It is not difficult to 
see that the orbit of a pair (x, y) under the action of G is characterized by the triple 

(|x|,|y|,|xny|). 

Following the notations of lITTl . the g-binomial coefficient [^] expresses the 
cardinality of X^. We have 

if g = 1 
if g > 1 



/ n-l 



n 
w 



n 

i=0 
n-l 

n 

i=0 



n 



w 
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In terms of the variable 



X 



if g = 1 



1 



we have 



ui— 1 1- 
w{n-w) "TT 1^ ~ ^ 



n 

i=0 



w — I 



- 1 

w{n—w) 



if Q > 1 ' 



n ! 



t« ! n — w ! 



We have the obvious decomposition into pairwise orthogonal G-invariant sub- 
spaces: 

C{X) = C{Xo) ± C(Xi) ±---± C{Xn). 
The decomposition of C(Xu,) into G-irreducible subspaces is described in ifTTl . 
We have 

where the Hf^^^ are pairwise isomorphic for equal k and different w. and pairwise 
non isomorphic for different k. The picture looks like: 



CiX)= CiXo)± C(Xi)± ... ±C(Xl|j)± 



^0,1 ^ 



±C(X„_i) ±C{Xn) 

-L Hi n-1 



where the columns represent the decomposition of C(Xt„) and the rows the isotypic 
components of C{X), i.e. the subspaces := ^ ± -fffc.fc+i -L • • • -L Hj^ n_]^, 
<k < [|J , with multiphcity rrik = {n - 2k + 1). 
Let, for all (A;, with 0<k<i<n — k, 

/ ^ V'fe,j(/) : ^fc,i(/)(2/) = E|x|=fc/(2;) 



and 



5k : CiXk) ^ C(Xfc_i) 

/ ^ (^fc(/):(^fc(/)(^) = EN=fc/(^) 



= ker 5k and 


= v^ 






k). Moreover, 


/ifc := dim(iJfc,fc) 




n 
k 




n 
k-l 





We need later the following properties of tljk,i- 
Lemma 5.2. Iff,gG Hk,k< 

\n - 2k' 



(16) 
Moreover, 
(17) 



l-k _ 

3 - k' 

i — k 



Mi-k) 



if, 9)- 
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Proof. ifTTl Theorem 3] proves (fT6l ). The relation (fTTl ) is straightforward: if \z 

^^,MkAm^) = E V'A.,(/)(2/) = E ( E /(^)) 

|3/|=i |s/|=« |a::|=A: 

|a;|=fc |j/|=* |a;|=A; 
xC2 xGyGz xGz 

V'fcj(/)(2)- 

□ 

Now we want to calculate the matrices of size = (n — 2A; + 1) associated 
to each isotypic space I^. We fix an orthonormal basis {ek^k,i, • • • ; ^k,k,h^.) of H^ k 
and we define '>Pk,i{^k,k,s)- It is clear from the definitions above that 

Cfc^i s can be assumed to take real values. From ([T6l ). for fixed A; and i, they form 
an orthogonal basis of H^ i with square norm equal to ["if ] g^^'"''). Normalizing 
them would conjugate Ek by a diagonal matrix, so we can omit to do it. The matrix 
Ek is indexed with i,j subject to A; < i,j < n — k. From the construction, we 
have Ek^ij{x,y) = if / i or \y\ / j; since the matrix E^ is zonal, we can 
define Pk,i,j by 

Ek,i,jix,y) = Pk,i,j{i - \xr\y\) 
and our goal is to calculate the Pk,i,j- It turns out that these functions express in 
terms of the so-called g-Hahn polynomials. 

We define the g-Hahn polynomials associated to the parameters n,i,j with < 
i < j < ntoht the polynomials Qk{n,i,j;x) with < k < min(i,n — j) 
uniquely determined by the properties: 

• Qk has degree k in the variable [x]. 

• (Qfe)fc is a sequence of polynomials orthogonal for the weights 

<u<i w{n,i,j;u) = 

• Qfc(o) = 1 

The polynomials defined in IITtII and 15.3.1] correspond up to multiplication by 
/ifc to the parameters {n,w,w) and, with the notations of ||T9l, according to The- 
orem 2.5, again up to a multiplicative factor, Qk{n,i, j; x) = Em{i,n — i,j,i — 
x; q~^). The combinatorial meaning of the above weights is the following: 

Lemma 5.3. |fT9l Proposition 3.1] Given x € Xi, the number of elements y € Xj 
such that |x n y| = i — u is equal to w{n, u). 

Theorem 5.4. If k < i < j < n — k, \x\ = i, \y\ = j, 

Ek,ij{x,y) = iXlhk^^^^j;^^ \xny\) 
If\x\ / ior\y\ / j, Ek^ij{x,y) = 0. 

Proof. We proceed in two steps: the first step (fTSl) calculates Pk,i,j{0) and the 
second step ([T9l ) obtains the orthogonality relations. 



fix) 



j - k 
i — k 



n — I 
j — i + u 



u{j—i+u) 
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Lemma 5.5. With the above notations, 



rj-fc"| rn-2fc 

(18) Pk,,m = \x\hu ^'~i\ 

Proof. We have Pk,i,j{^) = Ek^ij{x,y) for all x,y with |x| = i, \y\ = j, x C y. 
Hence 



kl L j-k 1 k(j-k) 



Pk,i,j{0) = Ek,ij{x,y) 



n R 

LjJ LiJ \x\=i,\y\=j 
xCy 



r P 

Iji 111 \x\=i,\y\=j s=l 
xCy 



r P 

LjJ 111 S = l \y\=j \x\=i 

xCy 

hk 



rmYUZ '^idi^k,i,s)iy)^k,j,s{y) 



[j] ih s=l \y\=j 



Since, from ([TtI) 



i^i,j{ek,i,s) = i^i.j ° i'k,i(ek,k,s 



j - k' 
i — k 



'4^k,j{^k,k,s) 



^k,jyS: 



we obtain 



k,i,j 



hk 



jj LiJ S=l \y\=j 



j - k 

i — k 



ek,j,s{y)^k,jAy) 



[j-k-\ hk ^j-k] \n-2k 



-^^|X|(efcj-„efcj-,) = \X\hk 

Vj\ Vi\ s=l 



k\ L j-k J k(j-k) 



from (fT6l ). 

Lemma 5.6. Wi'f/t above notations, 

i |^n-2fcj |-n-2fcj^fc(i+j-2fc) 

(19) Yw{n,i,j;u)Pk,iAu)Pl,i,j{'^)=h,l\X\^hk^—— ^~ 



□ 



u=0 



[•] 



Proof. We compute S := Eyex Ek,i,j{x,y)Ei^i'jf{y, 
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hk hi 

s = X] (^k,i,s{^)^k,j,s{y)^i,i' Ay)^i,j' Az) 

yeX s=l t=l 
hk hi 

s=l t=l yeX 
hk hi 



s=l t=i 
hk hi 



'^'^ek,i,s{x)eij>^t{z)\X\ 



=1 t=i 



^k,l^j,i'\X\ 



n - 2k' 

j - k 

n - 2U 
j - k _ 



'n - 2k 
_j-k 

hk 



qHj '')Y,ek,iAx)ei,f,s{z) 

s=l 

q''^^-^^Ek,^,r{x,z). 



We obtain, with j = i', j' = i, x = z e Xi, taking account of Eij i{y,x) 
Ei,i,jix,y), 



X] Ek,ijix,y)Ei^ij{x,y) = 6k,i\X\ 

y&Xj 



n - 2k' 
j - k _ 



q''^'-'^Ek,i,i{ 



x,x}. 



The above identity becomes in terms of Pj, 



k,i,j 



XI ^k,i,j{i' -\xn y\)Pi,ij{i -\xny\) = 6k,i\X 

yeXj 

Taking account of (ITSl ) and Lemma 1531 we obtain ( fT9l ) 



n - 2k' 
j - k _ 



I'^'-'^PkMo)- 



□ 



To finish the proof of Proposition 15 .41 it remains to prove that Pk,i,j is a polyno- 
mial of degree at most k in the variable [u] = [\xr\y\]. It follows from the reasons 
invoked in lITTl in the case i = j (see the proof of Theorem 5). □ 

Remark 5.7. In the case q = 1, i.e. the Hamming space, we could have followed 
the same line as for the sphere in order to decompose C{Hn) under the action of 
G. We could have started from the decomposition ofC{Hn) Q under the action 
ofT := T yi Sn = Aut(i/„) and then we could have decomposed each space 
Pk under the action of G = Stab(0",r). But we have a G-isomorphism from 
C(X^) = C( J-) to P^ given by: 

wt{y)=w 

Note that the inverse isomorphism is the Fourier transform on (Z/2Z)". So we 
pass from one to the other decomposition ofC{Hn) through Fourier transform. 
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6. An SDP upper bound for codes from positive definite functions 

In this section we want to explain how the computation of the continuous G- 
invariant positive definite functions on X can be used for applications to coding 
theory. In coding theory, it is of great importance to estimate the maximal number 
of elements of a finite subset C of a space X, where C is submitted to some con- 
straints. Typically X is a metric space with G-invariant distance d{x, y) and the 
constraints are related to the values taken by the distance on pairs of elements of 
C. In the following we concentrate on the basic case where the requirement is that 
the distance takes non zero values at least equal to some minimum 5. We denote 
by D the set of all values taken by d{x, y) and we define D>s = D r\[5, +oo[ and 

A{X, 5) := max{card(C) : d{c, c') > 6 for all c / c', (c, c') G C^}. 

We first focus on an upper bound for A{X, 6), which is obtained very obviously 
from the optimal value of the following program: 

Definition 6.1. 

m{X,6) =inf{ t: F e C{X^), F = F, F h 

(20) F{x,x)<t-1, 

Fix,y)<-1 d{x,y)>6} 

Then we obtain an upper bound for A{X, 6): 

Theorem 6.2. 

A{X,S) < m{X,5). 
Proof. For a feasible solution F, and for C C X with d{C) > 6 we have 

0< Yl F{c,c')<{t-l)\C\-\Cm-l) 
(c,c')eC2 

thus \c\<t. n 

Now the group G comes into play. From a feasible solution F one can construct 
a G-invariant feasible solution F' with the same objective value: 

F'{x,y) = / F{gx,gy)dg 
Jg 

thus we can add to the conditions defining the feasible solutions of m{X, 6) that F 
is G-invariant. Then we can apply Bochner characterization of the G-invariant pos- 
itive definite functions (Theorem 14.111 ). Moreover we have also seen in Theorem 
14.121 that if X is a homogeneous space, the finite sums of type ([T2l ) are arbitrary 
close for II II oo to the G-invariant positive definite functions on X, so we can replace 
F by an expression of the form ([T2l ) in the SDP m{X, 6). Moreover, we replace 
Ek{x, y) with its expression Yk{u{x, y)) in terms of the orbits of pairs and we take 
account of the fact that F = F. All together, with the notations of subsection 14.31 
we obtain the (finite) semidefinite programs: 

m^'^\X,5) =inf{ t: Fq ^ 0, . . . , Ffc ^ 0, . . . 

(21) Y,k>o{Fk,Yk{u{x,x))<t-l, 

Y.k>oiFk,Yk{u{x, y)) < -1 d{x, y) > 5] 
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where the matrices are real symmetric, with size nid^k, and Y}^{u{x,y)) = 
Yk{u{x, y))+Yk{u{x, y)). We insist that in the above program only a finite number 
of integers k are to be taken account of because ^ 7^ for a finite number of 
integers k. Thus we have 5) < m^'^\X, 6) and 

lim m^'^\x,5) = m{X,5). 

d— )-+oo 

6.1. The 2-point homogeneous spaces. We recall that a sequence of orthogonal 
functions {Pk)k>o is associated to X such that the G-invariant positive definite 
functions have the expressions 

F{x, y) = Y^ fkPk{d{x, y)) with h > 0. 

fe>0 

Then 

miX,6)=mf{l + Zk>ifk ■■ fk>0, 

1 + Ek>i fkPkii) < for all i£D>s} 

We restate Theorem 16.21 in the classical form of Delsarte linear programming 
bound: 

Theorem 6.3. LetF{t) = fo+fiPi{t) + - ■ ■ + fdPd{t)- Iffk > far all < k < d 

and /o > 0, and if F{t) < Ofar all t S Dyg, then 

A(X,<5)< + + + 
Jo 

Example: X = S''', d{x,y) = 6{x,y), d{C) = 7r/3. This value of the minimal 
angle corresponds to the kissing number problem. A very good kissing configura- 
tion is well known: it is the root system E^, also equal to the set of minimal vectors 
of the £'§ lattice. It has 240 elements and the inner products take the values ±1, 0, 
±1/2. We recall that the zonal polynomials associated to the unit sphere are pro- 
portional to the Gegenbauer polynomials in the variable x ■ y. \f P{t) obtains 
the tight bound 240 in Theorem [631 then we must have P(t) < for t G [-1, 1/2] 
and P{—1) = P(±l/2) = P(0) = (as part of the complementary slackness 
conditions). The simplest possibihty is P = (t - l/2)t^{t + l/2f'{t + 1). One 
can check that 

320 „ o 16 o 200 o 832 o 1216 o 5120 o 2560 o 
P = p8 H p8 H p8 H p| H p8 H p8 H p8 

3 ° 7 ^ 63 ^ 231 ^ 429 ^ 3003 ^ 4641 ^ 

and that 

^ = 240. 

/o 

Thus the kissing number in dimension 8 is equal to 240. This famous proof is 
due independently to Levenshtein |[27]| and Odlysko and Sloane |[35]| . A proof of 
uniqueness derives from the analysis of this bound (| 10]). For the kissing number 
problem, this miracle reproduces only for dimension 24 with the set of shortest 
vectors of the Leech lattice. For the other similar cases in 2-point homogeneous 
spaces we refer to ll28ll . 

It is not always possible to apply the above "guess of a good polynomial" 
method. In order to obtain a more systematic way to apply Theorem 16.31 one 
can of course restrict the degrees of the polynomials to some reasonable value, but 
needs also to overcome the problem that the conditions F{t) < for t G [—1, 1/2] 
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represent infinitely many linear- inequalities. One possibility is to sample the inter- 
val and then a posteriori study the extrema of the approximated optimal solution 
found by an algorithm that solves the linear program with finitely many unknowns 
and inequalities. It is the method adopted in [35], where upper bounds for the kiss- 
ing number in dimension n < 30 have been computed. We want to point out that 
polynomial optimization methods using SDP give another way to handle this prob- 
lem. A polynomial Q{t) G is said to be a sum of squares if Q = ^21=1 Qi 
some Qi G Being a sum of squares is a SDP condition since it amounts to 
ask that 



Here k is an upper bound for the degrees of the polynomials Qi. Now we can relax 
the condition that < for t G [-1,1/2] to F{t) = -Q{t) - Q'{t){t + l){t- 
l/2) with Q and Q' being sums of squares. A theorem of Putinar claims that in 
fact the two conditions are equivalent (but the degree of the polynomials under the 
squares are unknown). 

A very nice achievement of the linear programming method in 2-point homo- 
geneous spaces is the derivation of an asymptotic upper bound for the rate of 
codes (i.e. for the quotient logcard(C)/ dim(X)) obtained from the so-called 
Christoffel-Darboux kernels. This method was first discovered for the Hamming 
and Johnson spaces QUI and then generalized to the unit sphere |[24l and to all other 
2-point homogeneous spaces [28 1. It happens to be the best known upper bound 
for the asymptotic range. In [24] an asymptotic bound is derived for the density of 
sphere packings in Euclidean space which is also the best known. 

6.2. Symmetric spaces. For these spaces, which are not 2-point homogeneous, 
there may be several distance functions of interest which are G-invariant. For 
example, the analysis of performance of codes in the Grassmann spaces for the 
MIMO channel |[T4l involves both the chordal distance: 



and the product pseudo distance (it is not a distance in the metric sense): 



The reformulation of Theorem [6]2] leads to a theorem of the type [63] for any sym- 
metric function of the yi := cos 9i with the Jacobi polynomials Pfj.{yi, ■ ■ ■ ,ym) 
instead of the P^. For a general symmetric space, a theorem of the tvpe 16.3! is 
obtained, where the sequence of polynomials Pk{t) is replaced by a sequence of 
multivariate polynomials, and the set Ds is replaced by some compact subspace of 
the domain of the variables of the zonal functions, i.e. of the orbits of G acting 
on pairs. Then one can derive explicit upper bounds, see [45l for the permutation 
codes, \2i for the real Grassmann codes, [,37,1 and [141 for the complex Grassmann 
codes, |[T5l for the unitary codes, [|9l and fSTl for the ordered codes. Moreover an 
asymptotic bound is derived in [|2l and [|9l. 



Q = (1, t, . . . , t'')F{l, t,..., t'')* with FhO- 




m 



dp{p,q) ■■= Y[smei{p,q). 



i=l 
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6.3. Other spaces with true SDP bounds. An example where the bound (l20l) 
does not boil down to an LP is provided by the spaces V{n,q) endowed with the 
distance ([T5l) for which the matrices Ef^ are computed in section [533] (see [7]). In 
this case the group G is the largest group that acts on the SDP. 

Indeed, it is useless to restrict the symmetrization of the program (l20l ) to some 
subgroup of the largest group G that preserves {X, d). However, another interest- 
ing possibility is to change the restricted condition d{x, y) > S in A{X, 6) for the 
conditions: 

(22) d(x,y) > 6, d{x,e) < r, d{y,e) < r 

where e G X is a fixed point. Then the new A{X, e, r, 5) is the maximal number 
of elements of a code with minimal distance 6 in the ball B{e,r) C X. Here the 
group that leaves the program invariant is Stab(e, G). The corresponding bounds 
for codes in spherical caps where computed in f6l using the expressions of the 
zonal matrices of l5.5.1l 

We end this section with some comments on these SDP bounds. We have indeed 
generalized the framework of the classical LP bounds but the degree of understand- 
ing of the newly defined bounds is far from the one of the classical LP bounds after 
the work done since 1 17 1, see e.g. |28 1. It would be very interesting to have a better 
understanding of the best functions F that give the best bounds, to analyse explicit 
bounds and to analyse the asymptotic range, although partial results in these direc- 
tions have already been obtained. The fact that one has to deal with multivariate 
polynomials introduces great difficulties when one tries to follow the same lines as 
for the classical one variable cases. A typical example is provided by the config- 
uration of 183 points on the half sphere that seems numerically to be an optimal 
configuration for the one sided kissing number, and for which we failed to find the 
proper function F leading to a tight bound (see Uil). 

7. LOVASZ THETA 

In this section we want to establish a link between the program (l20l ) and the 
so-called Lovasz theta number. This number was introduced by Lovasz in the sem- 
inal paper |29| in order to compute the capacity of the pentagon. This remarkable 
result is the first of a long list of applications. This number is the optimal solution 
of a semidefinite program, thus is "easy to calculate", and offers an approxima- 
tion of invariants of graphs that are "hard to calculate". Since then many other 
SDP relaxations of hard problems have been proposed in graph theory and in other 
domains. 

7.1. Introduction to Lovasz theta number. A graph F = {V, E) is a finite set V 
of vertices together with a finite set E of edges, i.e. E C V"^. An independence 
set 5* is a subset of V such that 5^ H £^ = 0. The independence number aiV) 
is the maximum of the number of elements of an independence set. It is a hard 
problem to determine the independence number of a graph. The connection with 
coding theory is as follows: a code C of a finite space X with minimal distance 
d{C) > 5 is an independence set of the graph r(X, 5) which vertex set is equal to 
X and which edge set is equal to E^ := {{x,y) G X^ : d{x,y) G]0,d[}. Thus 
the determination of A{X, 6) is the same as the determination of the independence 
number of this graph. 
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Among the many definitions of Lovasz theta, we choose one which generalizes 
nicely to infinite graphs. For S C V, let be the characteristic function of S. Let 

M{x,y) ■.= ^^ls{x)ls{y). 

The following properties hold for M: 

(1) M G M"^", where \V\ = n, and M is symmetric 

(2) M ^ 

(3) E.evMix,x) = l 

(4) M{x,y) =0 if {x,y)£E 

(5) Ei.,y)eV^M{x,y) = \S\. 

Definition 7.1. The theta number of the graph T = [V, E) with y = {1, 2, . . . , n} 

is 

??(r) = max { Y.i j Bij : B G M"^", B^O 

(23) ' ' EiBi,i = l, 

Bij = {i,j)£E} 

The dual program for {} has the same optimal value and is equal to: 

T9(r) = min { t : BtO 

(24) Bi^i = t-l, 

Bi,j = -l {hj)iE] 

The complementary graph of T is denoted F. The chromatic number x(r) is the 
minimum number of colors needed to color the vertices so that no two connected 
vertices receive the same color. In other words it is a minimal partition of the vertex 
set with independence sets. Then the so-called Sandwich theorem holds: 

Tlieorem 7.2. 

«(r) < t9(f) < x(r) 

Proof. The discussion prior to the theorem proves the first inequality. For the sec- 
ond inequality, let c : y — )• { 1 , . . . , A;} be a coloring of F. Then the matrix C with 
Ci^j = — 1 if c{i) / c{j), Ci^i = k — 1 and Cjj = otherwise provides a feasible 
solution of □ 

7.2. Symmetrization and tlie (/-gones. Now we assume that G is (a subgroup of) 
the automorphism group Aut(F) of the graph. Then, G acts also on the above 
defined semidefinite programs. Averaging on G allows to construct a G-invariant 
optimal feasible solution B' from any optimal feasible solution B with the same 
objective value: 



\G\ 

Thus one can restrict in the above programs to the G-invariant matrices. Then 
one can exploit the method developed in previous sections, in order to obtain a 
description of the G-invariant B ^ form the decomposition of the space C{V) 
under the action of G. We illustrate the method in the case of the (/-gone Cq. There 
we have V = G = Zq the group of integers modulo q. Let be a fixed primitive 
root of 1 in C. Let Xk '■ "^g ^ C* be defined by Xkix) = Cg^- The characters of 
Zg are the Xk ^ ^ k < q — I and we have the decomposition 

c{Zq) = elilcxk- 
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According to Theorem 14.111 the G-invariant positive definite functions on V are 
exactly the functions F{x, y) of the form: 

9-1 q-l 

Fix,y) = Y.hXk{x)xk{y) = J^/fcC?"-^) 

fc=0 fc=0 

with fk>0- The ones taking real values have the form 

L9/2J 

F {x, y) = fk cos{{x -y)2k7r/q), fk>0. 

k=0 

When one replaces in the expression Bi j = F{i,j), the SDP transforms into a 
LP on the variables /fc. More precisely, we compute J2{xy)£V^ P{x,y) = q^fo 
and Ylxev -^(^' x) = Q Ylk fk- Thus we obtain (after a change of qfk to fk)'- 

^{Cg) =max{qfo: fk>0,0<k< [q/2\ , 

Vim 

k=0 

La/2J 

E /fcCOs(2/c7r/g) = 
fe=o 

The optimal value of this very simple linear program, is obtained for /i = /2 = 
• • • = fVq/2\-i = 0, and equals 



'^iS-'o) ~ \ qcos{n/q) 
I l+cos(7r/(jr) 



if q is even 
if q is odd . 



Note that when q is even, the independence number of the g-gone is exactly q/2. If 
the independence number of a graph as simple as the g-gone is not a great deal (it 
is of course equal to [q/2\), a more challenging issue is to determine its capacity. 
In general, the capacity C{T) of a graph T is defined to be 



C(r) = lim a(r 



.1/n 



-oo 



Here the graph is defined as follows: its vertex set is equal to V"- and an edge 
connects (xi, . . . , x„) and (yi, . . . , y„) iff for all 1 < i < n either Xi = yi or 
{xi, yi) ^ E. Introduced by Shannon in 1956, this number represents the effective 
size of an alphabet used to transmit information through the channel associated 
to the graph F (where two symbols are undistinguable if they are connected by 
an edge). If the capacity of a graph is in general very difficult to calculate, the 
theta number of a graph provides an upper bound for it because ??(F") = ■diV)'^ 
(see 1291 ). This upper bound is an equality for the pentagon since on one hand 
^(Cs) = \/5 from our previous computation, and on the other hand it is easy to 
see that a{{C^Y) = 5 (while a(C5) = 2); this is the way taken by Lovasz in [i29ll 
to prove that C{C^) = \/5. The determination of the capacity of the g-gone for q 
odd and greater than 5 is still opened. 
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7.3. Relation with Delsarte bound and with m{X, 5). We introduce a slightly 
stronger bound for a{T) with -d' and its dual form: 

^'{T) = max { Y.i j Bi,j : B^O, B>0 

(25) ' EiBi.i = l, 

Bi,,=0 {i,j)eE} 

-d'iV) = min { t: BtO 

(26) Bi^i < t-l, 

Bi,j < -1 (ij) i E] 

Since M(x^ y) > 0, we still have that a{T) < 'i?'(r). Again one can restrict in the 
above programs to the G-invariant matrices. It was recognized independently by 
McEliece, Rodemich, Rumsey, and Schrijver 139] that Delsarte bound of Theorem 
16.31 for A{Hn, 6) is equal to ??' for the graph T{X, 6), once the feasible set is re- 
stricted to the Aut(i7„)-invariant matrices, and similarly for the other finite 2-point 
homogeneous spaces. Indeed, by virtue of Theorem 14.111 the matrices B turn to 
be of the form B{x,y) = Ylik>o fkPk{d{x,y))- This symmetrization process is 
of great importance, not only because it has the great advantage to change an SDP 
to an LP, but also because it does change the complexity of the problem. Indeed, 
there are algorithms with polynomial complexity that do compute approximations 
of the optimal value of SDP's, thus algorithms with polynomial complexity in the 
number of vertices of F for d. But the graphs arising from coding theory have in 
general an exponential number of vertices, e.g. 2" for the Hamming graph. It is 
important to insist that the symmetrized theta has polynomial complexity in n. 
Now we can see that the program m{X, 6) (l20l ) is a natural generalization of 
for metric spaces under the assumptions of Section |4l We refer to [8] for a 
more general discussion about generalized theta where also chromatic numbers are 
involved. 

8. Strengthening the LP bound for binary codes 

In this section we explain how the zonal matrices Ek{x, y) related to the binary 
Hamming space computed in l5.5.3l are exploited in [40 1 in order to strengthen the 
LP bound. We shall work with the primal programs so we start to recall the primal 
version of (l20l ) in the case of the Hamming space. 

We recall that the sequence of orthogonal functions {Pk)o<k<n with = Kj. 
the Krawtchouk polynomials is associated to if„ such that Pi^{d{x, y)) ^ 0. As a 
consequence, we have for all > 

Pk{d{c,c'))>0. 

We introduce the variables Xj, for i G [0 . . . n] 

(27) Xi := -card{(c, c) € : d{c,c') = i}. 

card(C j 

They satisfy the properties: 

(1) xo = 1 

(2) Xi>0 

(3) EiXiPkii) > OforalU > 

(4) Xi = if i e [1 . . . 5 - 1] 
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(5) card(C) = ^ • Xj. 
With these properties which are linear inequalities, we obtain the following linear 
program which is indeed the dual of (l20l ): 

sup { 1 + Yh=5 ■■ Xi> 0, 

1 + E»=5 XiPkii) > for all 1 < A; < n } 

where we have taken into account Pq = 1. 

We recall that to every < k < [|J , we have associated a matrix Ei.{x, y) h 
of size n — 2fc + 1. In particular, for all C C Hn (see the remark l4~T0l) . 

Ek{c,c')hO. 

These constraints are not interesting for pairs because they are not stronger than 
the linear inequalities coming from the Krawtchouk polynomials. They are only 
interesting if triples of points are involved: namely we associate to (x, y, z) G H"^ 
the matrices 

Fk{x,y,z) := Ek{x - z,y - z). 
We have for all C C Hn, and for all z G Hn, 

J2 Fk{c,c',z)hO 

which leads to the two positive semidefinite conditions: 

Eic,c',c")ec^Pk(.c,c',c")hO 
E{c,c')ec2, c"<^cPk{c,c',c") t 

Theorem 15 .41 expresses the coefficients of Ek{x — z,y — z) in terms of wt{x — z), 
wt{y — z), wt{x — y); so with a := d{y, z), b := d{x, z), c := d{x, y), we have 
for some matrices Tk{a, b, c), 

Fk{x,y,z) = Tk{a,b,c). 

We introduce the unknowns Xa,b,c of the SDP. Let 

a + 6 + c = mod 2 
a + b + c<2n 
n:= ^{a,b,c) e[0...nf : c<a + b | 

b < a + c 
a <b + c 

It is easy to check that = {{d{y, z),d{x, z),d{x, y)) : {x, y, z) G H^}. Let, for 
(a, 6, c) G Q., 

Xa,b,c ■■= JTT^T card{(x, y, z) e : d{y, z) = a, d{x, z) = b, d{x, y) = c}. 

card(G) 

Note that 

xo,c,c = -77;rcard{(x,y) G : d{x,y) = c} 

card(C j 

thus the old variables Xi (|27] ) of the linear program are part of these new variables. 
We need a last notation: let 

i(a, 6, c) := card{z G Hn ■ d{x, z) = b and d{y, z) = a} for d(x, y) = c 
= (?)(n) wherea-& + c = 2i 



(28) 
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Then, if C is a binary code with minimal distance at least equal to 6, the following 
inequalities hold for Xa,b,c '■ 

(1) a;o,o,o = 1 

(2) Xa,b,c > 

(3) Xa,b,c = Xr{a),T{b),T{c) for all permutation r of {a, b, c} 

(4) Xa,b,c < t{a, b, c)xo,c,c, Xa,b,c < t{b, C, a)xo^a,a, Xa,b,c < t{c, a, 6)xo,6,fe. 

(5) Ea,b.c Tk{a, b, c)xa,b,c ^ for all < fc < [fj 

(6) Ea,b,c Tk{a, b, c){t{a, b, c)xo,c,c - Xa,b,c) ^ for all < A; < [f J 

(7) Xa^b,c = if a, 5 or c g]0, 6[. 

(8) card(C) = ^o,c,c- 

Conditions (5) and (6) are equivalent to (l28l) . Condition (7) translates the assump- 
tion that d{C) > 6. Thus an upper bound on card(C) is obtained with the optimal 
value of the program that maximizes X^f,xo,c,c under the constraints (1) to (7). 
This upper bound is at least as good as the LP bound because the SDP program 
does contain the LP program of 16.11 Indeed, the sum of the two SDP conditions 
(|28] ) is equivalent to 

^ Ek{x - z,y - z) hO- 

zeHn 

We claim that this set of conditions when k = 0, 1, ... , [^J is equivalent to the 
set of conditions Pk{d{x,y)) >z for A; = 0, . . . , n. Indeed let Bj^{x,y) := 
Y.z&Hr, ^k{x - z,y - z). Up to a change of Bk{x,y) to ABk{x,y)A*, we as- 
sume that Ek was constructed using the decomposition of C{Hn) first under F := 
T X S'n = Aut(i^„) then under G (see Remark 157/] ). Clearly Bk is F-invariant. 
Since x — )• Ek^ij{x,y) S Pi and Pi is a F-module, also x — )• Bk,ij{x,y) S Pi 
and similarly y — )• Bf^^ij{x,y) G Pj. But Pi and Pj are non isomorphic F- 
modules for i j thus Bk^ij{x,y) = for i / j. Since Pi is F-irreducible, 
Bk,i,i{x, y) = XiPi{d{x, y)) for some Aj > that can be computed with Bk{x, x). 
So we have proved that the linear program associated to Hn like in l6.1l is contained 
in the SDP program obtained from the above conditions (1) to (7). Moreover it 
turns out that in some explicit cases of small dimension the SDP bound is strictly 
better than the LP bound (see BOl ). 

A similar strengthening of the LP bound for the Johnson space and for the spaces 
of non binary codes where obtained in 140] and [20,1 . In the case of the spherical 
codes, for the same reasons as for the LP bound, one has to deal with the dual 
program, see fill . 
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